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For : Method and Apparatus for Measuring a Device Under Test 
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Petition For Acceptance of Specification filed in 
Response to Notice to File Corrected Application 

Papers 

Dear Sir, 

Applicant petitions under 37 CFR §1.181 and 
MPEP §1002.02 (b) (3) to invoke the supervisory 
authority of the Commissioner in matters not 
otherwise expressly provided for. 



Facts Upon Which Petition is Presented 
Petitioner filed the present patent application 
on August 5 th , 2003. On November 3 rd , 2003, The USPTO 
mailed a Notice to File Corrected Application Papers 
indicating that the Specification pages 43-93 (the 
last 50 pages) were not in compliance with 37 CFR 
§1.52 (a). See Exhibit A. The cited text of 37 CFR 
§1.52 (a) found in the Notice reads: 

Papers must be legibly written either by a 
typewriter or a mechanical printer in permanent 
ink or its equivalent in portrait orientation 
on flexible, strong, smooth, nonshiny, durable, 
and white paper. Application papers must be 
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presented in a form having sufficient clarity 
and contrast between the paper and the writing 
thereon to permit the direct reproduction of 
readily legible copies in any number by use of 
photographic, electrostatic, photo-offset, and 
microfilming processes and electronic 
reproduction by use of digital imaging and 
optical character recognition. 

The copy of the Specification held by 
Undersigned appears to be in full conformance with 
the cited rule. Accordingly, Undersigned telephoned 
the Customer Service Center of the Initial Patent 
Examination Division to obtain additional specifics 
as to what was not in conformance in order to 
properly respond to the Notice. Undersigned was 
told that, although, additional specifics were not 
available, the cause of the Notice was most likely a 
problem with the weight of the paper upon which the 
Specification was printed. Please note that a 
review of the actual filed specification was not 
performed by the USPTO representative. 

In reliance on this information, on November 
17 th , 2003, Petitioner filed a Response to the Notice 
together with a clean copy of the Specification 
using Hammermill bond, A4 paper. This paper is high 
quality and is the paper that has been used by 
Undersigned and others in filing Applications and 
Responses to the USPTO without incident for many 
years and is the same kind of paper used for the 
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present Petition. A copy of the second-filed 
Specification is attached hereto as Exhibit B. 

On December 22 nd , 2 003, the USPTO sent a second 
Notice to File Corrected Application Papers citing 
the same rule and pages. See Exhibit C. 
Undersigned again called the Customer Service Center 
of the Initial Patent Examination Division to obtain 
additional specifics as to what might be the 
problem. The USPTO Representative mentioned that it 
appeared on her computer display that a black line 
showed on the pages cited. Once again, no review of 
the actual filing was performed by the USPTO 
Representative. A review of copies of the 
Specification held by Undersigned and the Assignee 
failed to uncover the black bars claimed by the 
Initial Patent Examination Division. 

Specific Request by Petitioner 
Applicant petitions the Commissioner for 
acceptance of the second filed copy of the subject 
application. In the alternative, if the second copy 
as filed has been destroyed, the Application 
petitions the Commissioner for acceptance of the 
attached copy (Exhibit B) of the second-filed 
application as filed on November 3, 2003. 

If Commissioner determines that the second- 
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filed Specification does not conform to the cited 
rule, Petitioner further petitions to receive 
instructions of specific changes that must be made 
to the Specification to bring it into conformance 
with the cited rule . 

So as to avoid abandonment during the petition 
process, Applicant has responded to the outstanding 
Notice, issued December 22, 2003, with yet another 
copy of the Specification. 

If any clarifications can be made by way of 
telephonic interview, Commissioner's delegate is 
invited to contact the Undersigned. 



Respectfully submitted, 



Vane Adamian 
Petitioner 




/ F^tent Attorney 
/Registration No. 37, 928 
Attorney for Petitioner 
Phone: (970) 206-9177 
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Method and Apparatus for Measuring a Device Under Test 

Using an Improved 

Through-Ref lect-Line Measurement Calibration 

Background 

[01] U.S. Patent Application Ser. No. 10/294,317 filed 
November 14, 2002 entitled "Method and Apparatus for 
Performing Multiport Through-Ref lect-Line Calibration 
and Measurement" (herein "the '317 Patent Application") 
teaches a method and apparatus for performing a 
through-ref lect-line calibration with a multiple port 
vector network analyzer ("VNA"). The disclosed method 
teaches determination of ten-term error coefficients 
that may be used in a method disclosed in US Patent 
5,434,511 (herein "the '511 Patent") to arrive at S- 
parameters of a device under test (herein "the DUT") 
with systematic errors attributable to the VNA and the 
multiport test set corrected out of the measured S- 
parameters of the DUT. In the method according to the 
'317 Patent Application, it is necessary to know a type 
of high reflect calibration standard used in the 
process. Also, in the method disclosed in the '317 
Patent Application, if the through calibration standard 
used in the calibration procedure is a non-zero 
through, then its electrical length must be known and 
specified by the calibration standard manufacturer in 
order to fully characterize the error coefficients in 
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the calibration process. In many cases, the electrical 
length of the non-zero thru is either unknown, not 
known to a sufficient accuracy, or is a non-repeatable 
value. A through calibration standard is typically 
unknown for on-wafer calibrations because the landing 
of the probe to the calibration standard is non- 
repeatable, which renders the electrical length also 
non-repeatable, and therefore, unpredictable. Even 
though manufacturers of coaxial calibration standards 
often include electrical length data, over time and 
after multiple uses, the electrical length can change 
making a calibration process that relies on the 
electrical length specifications, less accurate. Due 
to machining, process, and other non-repeatable 
manufacturer variations, the manufacturer 
specifications of the through calibration standards may 
not be as accurate as the calibration and measurement 
presumes. In addition, frequent use of the through 
calibration standard can cause the dimensions and, 
therefore, the calibration specifications to change 
over time while the calibration process presumes 
otherwise. In high frequency on-wafer measurements, a 
non-zero through is used exclusively. Furthermore, in 
an on-wafer calibration process, it is difficult to 
land probes repeatably on a calibration substrate to 
coincide with the electrical delay value provided by 
the calibration standard manufacturer. Any discrepancy 
between the actual electrical delay and the presumed 
electrical delay causes uncertainties in the 
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calibration and, therefore, the corrections made to the 
measurement of a device under test ( "DUT" ) The 
uncertainties are exacerbated as measurement 
frequencies increase . 

[02] There is a need, therefore, to improve the through- 
reflect-line calibration process to accommodate a thru 
calibration standard with an unknown delay or an 
unknown type of high reflect calibration standard with 
a known phase offset. 



Brief Description of the Drawings 

[03] Figure 1 is a diagram of a vector network analyzer 
("VNA"), switch network and signal separation devices 
used for measurement of a device under test ("DUT"). 

[04] Figure 2 shows the VNA with a reflect calibration 
standard connected . 

[05] Figures 3 and 4 show the VNA with a line calibration 
standard connected and respective switch network 
configurations for forward and reverse measurements. 

[06] Figure 5 shows the VNA with a matched load 
calibration standard connected. 

[07] Figures 6 and 7 show the VNA with a source 

terminated through calibration standard connected and 
respective switch network configurations for forward 
and reverse measurements. 
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[08] Figures 8 and 9 show the VNA with a locally 

terminated through calibration standard connected and 
respective switch network configurations for forward 
and reverse measurements . 

[09] Figure 10 is a flow diagram of some of the error 
coefficients determined in a method according to the 
present teachings . 

[10] Figures 11 and 12 show an optional VNA measurement 
with a locally terminated through calibration standard 
connected between one of the proximal pairs of 
measurement ports and respective switch network 
configurations for forward and reverse measurements. 

[11] Figures 13 through 17 is a flow chart of a method 
for calibrating and correcting S-parameters of a 
measured DUT according to the present teachings. 

Detailed Description 

[12] US Patent No. 5,434,511 entitled "Electronic 
Microwave Calibration Device" (herein "the f 511 
Patent") teaches a method for correcting measured S- 
parameters to arrive at actual S-parameters for a 
measured device. The method uses twelve error terms 
calculated from a short-open-load ("SOLT") calibration 
process. As one of ordinary skill in the art 
appreciates, a through-ref lect-line ("TRL") calibration 
process is also able to determine the twelve error 
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terms from which actual S-parameters are determined 
based upon measured S-parameters of a device under 
test. The twelve error terms are: (1) forward 
directivity ("EDF"), (2) forward source match ("ESF"), 
(3) forward reflection tracking ("ERF"), (4) forward 
isolation ( "EXR" ) , (5) forward load match ("ELF"), (6) 
forward transmission tracking ("ETF"), (7) reverse 
directivity ("EDR"), (8) reverse source match ("ESR"), 
(9) reverse reflection tracking ("ERR"), (10) reverse 
isolation ("EXR"), (11) reverse load match ("ELR"), and 
(12) reverse transmission tracking ("ETR") . Because 
the isolation of the switches used in a multiport test 
set is so high, it is more accurate in most 
circumstances to ignore the forward isolation ("EXF" ) 
and reverse isolation ("EXR") terms and set them to 
zero. This yields a 10-term error correction process. 
The TRL calibration process is more time consuming, but 
has the capability for more accurate measurements. 
More accurate measurements yield improved error terms 
and improved corrected S-parameters. The following 
discussion is directed toward multiple embodiments of 
an improved TRL calibration, DUT measurement, and S- 
parameter correction process. 

[13] With specific reference to FIGURE 1 of the 

drawings, there is shown a system block diagram of a 2- 
port VNA 200 in conjunction with a 4-port test set. 
The illustrated 4-port test set comprises a switch 
network 150 and two external signal separation devices 
114 2 and 114 4 , which may be used with a 2-port VNA 200 
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for measuring a 4-port device under test 101 (herein 
"the DUT") . Alternative embodiments of test sets 150 
having additional switches and signal separation 
devices may be used for measuring DUTs having more than 
4 ports. For purposes of the present disclosure , r the 
term "vector network analyzer (VNA) 200" refers to 
either a two port VNA alone, a 2-port VNA with a 
multiport test set for calibrating and measuring 
devices with more than 2 ports, or a multiport VNA. 
Before the DUT 101 is connected, a TRL calibration 
process is performed to determine systematic errors to 
a measurement as a result of the VNA 200 and associated 
connections. The VNA 200 deploys first and second 
reference channels 201, 202, respectively, and first 
and second reference channel samplers 110, 210, 
respectively. In the test set-up shown in FIGURE l,the 
samplers 110, 210 in a specific example may be bridges 
or directional couplers. The reference channel 
samplers 110, 210 are placed in the first and second 
signal paths 107, 108 on an opposite side of a signal 
transfer switch 106 from signal generator 105. The 
samplers 110, 210 extract in one direction a small and 
predictable portion of the signal present on the first 
and second signal paths 107, 108 for measurement by the 
first and second reference channels 201, 202, 
respectively. The sampled portion is typically -lOdB 
to -20dB for a coupler based sampler and approximately 
-6dB for a bridge based sampler relative to the signal 
level on the signal path 107 or 108. The source 



6 



Agilent PDNO: 10021239 
inventor: Vahe Adamian 



transfer switch 106 connects either the signal 
generator 105 to the first signal path 107 and a signal 
transfer switch terminating load 109 to the second 
signal path 108 or connects the signal generator 105 to 
the second signal path 108 and connects the signal 
transfer switch terminating load 109 to the first 
signal path 107. In a specific embodiment, there are 
only two switch positions for the signal transfer 
switch 106. 

[14] The illustrative example shown in the drawings and 
described herein includes a 4-port DUT 101 connected to 
measurement ports 103i, 103 2 , 103 3/ and 103 4 . The 
teachings herein disclose process details for 
calibrating ports 103i and 103 3 . The disclosed process 
is repeated for all remaining ports. The present 
teachings may also be applied to a multiport test set- 
up for measuring DUTs having two or four or more device 
connections. A switch network 150 permits connection 
of each measurement port 103i through 103 2 n to a first 
or second signal path 107, 108, the transfer switch 
terminating load 109, or to a local terminating 
impedance, 104i through 104 2N - Certain switch network 
configurations connect one of the measurement ports 103 
to the first signal path 107 or/and a different one of 
the measurement ports 103 to the second signal path 108 
while the remaining paths are terminated in the local 
terminating impedance 104 . The switch network 150 also 
has sampling arms 113, sampling arms 113i through 113 4 
in the illustrated embodiment. The sampling arms 113i 
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through 113 4 each comprise a sampler 114 that samples a 
small and predictable portion of the signal level 
present at the respective measurement port 103. The 
sampler 113 may be a coupler or a bridge that takes 
somewhere bet ween — 6dB and — 20dB of the signal level 
from the signal level present on the respective 
measurement port 103. In a specific embodiment 
according to the present teachings, the portion sampled 
from the measurement ports 103 is substantially the 
same portion sampled from the signal paths 107, 108. 
The sampled signal may then be connected to either a 
first or second test channel receiver 111, 112 through 
a respective sampling switch 115 or may be connected to 
a sampling arm terminating load 116. A switch network 
150 of this configuration may connect a reflection path 
from the measurement ports 103 to one of first and 
second test channel receivers 111, 112 while 
terminating the reflection paths from measurement ports 
103 not connected to a test channel receiver in a local 
sampling arm terminating impedance 116. 

[15] In an embodiment of a method according to the 
present teachings, a TRL calibration on a multi-port 
DUT having N device connections is performed by 
conducting a 2-port TRL calibration first on each one 
of N/2 direct pairs of the measurement ports 103 and if 
there are more than two measurement ports, on (N(N- 
2))/4 indirect pairs of the measurement ports 103. A 
user may define the N/2 direct pairs by representing 
all measurement ports 103 in groupings of two 
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measurement ports 103, where a first port in the direct 
pair is capable of connection to the first test channel 
receiver 111 and a second port in the direct pair is 
capable of connection to the second test channel 
receiver 112. As an example, if there are N 
measurement ports, the direct pairs of measurement 
ports are measurement ports 103i and 103 N/2 +i, 
measurement ports 103 2 and 103 N / 2+2 , et seq. up to the 
direct pair of measurement ports 103 N/2 and 103 N where 
measurement ports 103i through 103 N/2 are capable of 
connection to the first test channel receiver 111 and 
measurement ports 103( N+ i)/2 through 103 N are capable of 
connection to the second test channel receiver 112. 

[16] There are (N(N-2))/4 indirect pairs, which are 
defined as those unique groupings of two measurement 
ports 103 that are not represented in the set of direct 
pairs of measurement ports and where a first 
measurement port in the indirect pair is capable of 
connection to the first test channel receiver 111 and 
the second measurement port in the indirect pair is 
capable of connection to the second test channel 
receiver 112. In the illustrated example, there are 
two direct pairs: first direct pair comprising 
measurement ports 103i and 103 3 and second direct pair 
comprising measurement ports 103 2 and 103 4 . Also in 
the specific example, there are two indirect pairs: 
first indirect pair comprising measurement ports 103i 
and 103 4 and second indirect pair comprising 
measurement ports 103 2 and 103 3 . Measurement ports 103 
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capable of connection to the same test channel 
receiver, either the first test channel receiver 111 or 
the second test channel receiver 112, are herein 
referred to as "proximal pairs". An N-port DUT has 
N(N-2)/4 unique proximal pairs of measurement ports 
103. In the specific illustrated embodiment having 
four measurement ports, the proximal pairs are 
measurement ports 103i and 103 2 , measurement ports 103 3 
and 103 4 . 

[17] In the TRL calibration process and with specific 
reference to FIGURE 2 of the drawings, a high reflect 
calibration standard 301 ("the reflect 301") is 
measured at all measurement ports 103. The reflect 301 
may be either a short or an open circuit with or 
without an offset. In one embodiment of the process, 
it is not necessary to know the type of reflect 
standard used, but it is important to know the phase 
offset of the reflect 301 in order to calculate an 
electrical length of a thru calibration standard. In 
an alternative embodiment, it is not necessary to know 
the type of reflect or the phase offset of the reflect 
calibration standard as long as the electrical length 
of the thru calibration standard is known. In the 
illustrated example, the reflect 301 is connected to 
measurement port 103i. The switching network 150 is 
set so the measurement port 103i is connected to the 
first signal path 107 and the respective sampling arm 
113i is connected to the first test channel receiver 
111. All remaining unused measurement ports 103 2 , 103 3 
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and 103 4 are terminated in their respective local 
terminating loads 104 and their respective sampling 
arms 113 are connected to the sampling arm terminating 
loads 116 2 , II63, and 116 4 . As one of ordinary skill in 
the art appreciates, for measurement of the measurement 
port 103i only the switch network configuration that 
terminates in a characteristic impedance for those 
measurement ports 103 that are capable of connection to 
the first test channel receiver 111 are important to 
the results. Because the isolation of the switches 
that comprise the switch network 150 is so high, the 
measurement ports 103 capable of connection to the 
second test channel receiver 112 do not figure in the 
high reflect measurement that uses the first test 
channel receiver 111. Similarly, the measurement ports 
103 capable of connection to the first test channel 
receiver 111 do not figure in the high reflect 
measurement that uses the second test channel receiver 
112. The signal generator 105 is then swept through a 
desired frequency range as programmed by an operator 
and measurements are taken at specific frequencies 
across the range. During the frequency sweep the VNA 
200 measures a reference signal as well as a reflected 
signal and stores a ratio of the measured reflected 
signal level at the test channel receiver, 111 in the 
example, over the measured signal level at the 
reference receiver, 201 in the example. The system 
stores a resulting ratio, which is a frequency 
dependant reflection coefficient, herein referred to as 
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a high reflect characteristic for measurement port 
103a. 

A r eflect_l/Rlref lect_l 

[18] The same process of measuring the reflect 301 is 
repeated for each measurement port 103. For each 
measurement port, the switch network 150 is configured 
so that the measurement port 103 with the reflect 301 
receives the stimulus signal 105 and the switch network 
150 has a sampling arm 113 arrangement whereby the 
reflected signal is measured by the appropriate test 
port 111 or 112. The unused measurement ports 103 are 
terminated in the local characteristic impedance 104 in 
the switching network 150. The sampling arms 113 of 
the unused measurement ports 103 are also terminated in 
the respective sampling arm terminating loads 116. The 
signal generator 105 stimulates the appropriate signal 
path 107 or 108 with a signal that is swept over the 
same desired frequency range as in the reflection 
measurement of measurement port 103i. The VNA 200 
measures and stores a measurement of a ratio of the 
measured stimulus signal level over measured reflection 
signal yielding a frequency dependent array of 
reflected signal level over the stimulus signal, either 

A r eflect_n/Rlreflect_n Or B re f lectin / R2 re f lectin / Where "n M refers 

to the measurement port 103. 

[19] With specific reference to FIGURE 3 of the 
drawings, a next step in the calibration process 
connects and measures a low-loss delay line calibration 
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standard 401 ("the line 401") between each direct pair 
and indirect pair of measurement ports 103. The line 
401 may be an airline i.e. a delay line having an air 
dielectric, which is commonly used in metrology 
laboratories. Alternatively and as typical for on- 
wafer measurements, a delay line may be used. The 
delay of the line 401 is unknown, but the physical 
dimensions of the line 401 relate to a valid range of 
calibration frequencies. Additional delay line 
calibration standards can be used to cover a broader 
frequency range as desired. The delay of the line 401 
is defined over a range of frequencies bounded by 
approximately more than 20 degrees phase shift at a 
lowest specified frequency and less than 160 degrees 
phase shift at a highest specified frequency. At 
frequencies around 500 MHz and below, coaxial airline 
dimensions become very large and not practical. To 
measure at the lower frequencies, and with specific 
reference to FIGURE 5 of the drawings, two high-quality 
matched loads 501 ("the matched loads 501") are 
connected to each measurement port 103 of the direct or 
indirect pair being measured. The matched loads 501 
are used for calibrating the VNA 200 in a frequency 
range down to the lowest VNA frequency. The resulting 
calibration values of the line 401 and the matched 
loads 501 are different, but the algorithmic 
formulations using the measured ratios are the same. 

[20] With specific reference to FIGURE 3 of the 

drawings in which the line 401 is measured between the 
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first direct pair comprising measurement ports 103i and 
103 3/ the transfer switch 106 is set so that the signal 
generator 105 stimulates the first signal path 107 and 
the second signal path 108 is terminated at the source 
109. The switch network 150 is configured so that the 
measurement port 103i accepts the stimulus signal from 
the first signal path 107 and the signal from the 
sampling arm 113i is presented to the first test port 
111. The switch network 150 is further configured so 
that the measurement port 103 3 is terminated in the 
transfer switch characteristic impedance 109 through 
second signal path 108 and a transmitted signal is 
presented to the second test port 112 through sampling 
arm 113 3 . The signal generator 105 sweeps the desired 
frequency range and the VNA 200 measures signal level 
from the first and second test channel receivers 111, 
112 and the first and second reference receivers 201, 
202 and stores the results in data arrays. For 
purposes of clarity and consistency only, when the 
signal generator 105 is connected to the first signal 
path 107, all resulting measurements are considered 
forward direction measurements. Accordingly, the 
measurements made of the line 401 in the forward 
direction are represented as the data arrays: 

AfjLine_13/ 
Bf_iine_13/ 
Rlf_line_13/ and 
R2 f line 13 • 



y 
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where each array comprises a series of measured points 
at specific frequencies along the desired frequency 
range . 

[21] With specific reference to FIGURE 4 of the 
drawings, the line 401 remains connected and the 
transfer switch 106 is reconfigured so that the signal 
generator 105 stimulates the second signal path 108 and 
the first signal path 107 is terminated in the transfer 
switch characteristic impedance 109 at the source. The 
configuration of the switch networks 150 is not changed 
from the forward direction measurements. The signal 
generator 105 again sweeps the desired frequency range 
and the VNA 200 measures signal level from the first 
and second test channel receivers 111, 112 and the 
first and second reference receivers 201, 202 and 
stores them in data arrays. For purposes of clarity and 
consistency only, when the signal generator 105 is 
connected to the second signal path 108, all resulting 
measurements are considered reverse direction 
measurements. Accordingly, the measurements made of 
the line 401 in the reverse direction are represented 
as the data arrays: 

A r _line_13, 
B r _line_13 / 
Rlr line 13/ and 
R2 r _line_13 - 

where each array comprises a series of measured points 
at specific frequencies along the desired frequency 
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range . 

[22] If a broader frequency range is necessary, the 
same measurement procedure is performed on the first 
direct pair, measurement ports 103i and 1033 in the 
specific embodiment, with a different airline covering 
a different frequency band. In addition, the matched 
loads 501 as shown in FIGURE 5 of the drawings may be 
used to simulate a high loss line having a perfect 
match in order to take measurements at lower 
frequencies than are practical with an airline 
calibration standard. Depending upon the matched load, 
the quality of the match at higher frequencies, and the 
desired frequency range, the matched loads may be used 
in place of the airline calibration standard. As new 
measurements are made for the different frequency 
ranges using the appropriate calibration standards, the 
results are stored in the forward and reverse direction 
arrays with each data point corresponding to a specific 
stimulus signal frequency. Accordingly, the 
calibration frequency band can extend over more 
frequencies than is possible with a single airline 
calibration standard . 

[23] With specific reference to FIGURE 6 of the 

drawings, a next step in the calibration process is to 
connect and measure a source terminated through 
calibration standard 601 ("the thru 601") between each 
direct pair and if there are more than two measurement 
ports 103, each indirect pair. The thru 601 may have 
either a zero length or a non-zero length and in one 
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embodiment, it is not necessary that the electrical 
length be a known value. In an alternate embodiment, 
if the electrical length of the thru 601 is known, it 
is not necessary to know the phase offset of the high 
reflect calibration standard 301. For on-wafer 
measurements, it is not possible to obtain a high 
quality zero thru calibration standard. Accordingly, 
for on-wafer measurements, a non-zero thru calibration 
standard and zero offset high reflect calibration 
standard is typically used. 

[24] The first measurement of the thru 601 is made of a 
source terminated thru 601. Using the first direct 
pair as an example, the transfer switch 106 is set so 
that the signal generator 105 stimulates the first 
signal path 107 and the second signal path 108 is 
terminated in the transfer switch characteristic 
impedance 109. The switch network 150 is configured so 
that the measurement port 103i accepts the stimulus 
signal from the first signal path 107 and the sampling 
arm 113i is connected to the first test port 111. The 
switch network 150 is further configured so that the 
measurement port 1033 is terminated in the transfer 
switch characteristic impedance 109 through second 
signal path 108 and the sampling arm 1133 is connected 
to the second test port 112. The unused measurement 
ports 103 are terminated in the local characteristic 
impedances 104 and the unused sampling arms 113 are 
terminated in local sampling arm terminating loads 116. 
The signal generator 105 sweeps the desired frequency 
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range and the VNA 200 measures the signal level from 
the first and second test channel receivers 111, 112 
and the first and second reference receivers 201, 202 
and stores the results in memory. According to the 
nomenclature used for purposes of the present 
disclosure, because the signal generator 105 is 
connected to the first signal path 107, the resulting 
measurements are considered forward direction 
measurements. Accordingly, the measurements made of 
the source terminated thru 601 in the forward direction 
are represented as the arrays: 

Afnm_thru/ 
B f nm thru/ 
Rlfnm_thru/ ^nd 
R2 f nm_thru • 

where each array comprises a series of measured points 
at specific frequencies along the desired frequency 
range and subscripts "n" and "m" refer to the 
measurement ports between which the thru 601 is 
connected where "m" is the source port and "n" is the 
receive port. 

[25] With specific reference to FIGURE 7 of the 
drawings, the thru 601 remains connected and the 
transfer switch 106 is set for the reverse 
measurements. Specifically, the signal generator 105 
stimulates the second signal path 108 and the first 
signal path is terminated in the transfer switch 
characteristic impedance 109. The switch network 150 
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is unchanged. The signal generator 105 again sweeps 
the desired frequency range and the VNA 2 00 measures a 
signal level from the first and second test channel 
receivers 111, 112 and the first and second reference 
receivers 201, 202 and stores them in memory. Because 
the signal generator 105 is connected to the second 
signal path 108 , the resulting measurements are 
considered reverse direction measurements. 
Accordingly, the measurements made of the source 
terminated thru 601 in the reverse direction are 
represented as the arrays : 

A rn m_thru/ 
B r nm_thru/ 
Rlrnm_thru/ and 
R2 rnm thru - 



where each array comprises a series of measured points 
at specific frequencies along the desired frequency 
range . 

[26] At this point in the process, enough measurements 
have been taken to calculate forward and reverse 
directivity, source match and reflection tracking for 
all direct and indirect pairs. The present 
organization of the process, however, completes all 
measurements and then calculates all ten error terms 
for each direct, indirect, and proximal pairs of 
measurement ports. It is further possible with the 
measurements taken thus far to calculate an electrical 
length of the sum of the thru calibration standard 601 
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and zero or offset high reflect standard 301 that is 
used in the calibration process by extrapolating an 
argument of the reflection tracking term to DC. It is 
also possible to determine a type of high reflect 
standard used in the calibration process. Additional 
details of the specific error term and electrical 
length calculations are described herein. 

[27] With specific reference to FIGURE 8 of the 

drawings and with the thru 601 still connected, the 
transfer switch 106 is configured so that the signal 
generator 105 is in the first signal path 107 and the 
second signal path 108 is terminated in the 
characteristic impedance 109. The switch network 150 
is configured so that the measurement port 103 is 
connected to the first signal path 107 and the 
respective sampling arm 113 is connected to the first 
test channel receiver 111. The remaining unused 
measurement ports 103 capable of connection to the 
first test channel receiver 111 are terminated at 
respective local characteristic impedances 104. In 
addition, the sampling arms 113 of the unused 
measurement ports 103 are terminated in local sampling 
arm characteristic impedances 116. The switch network 
150 is further configured so that the measurement port 
103 capable of connection to the second test channel 
receiver 112 in the first direct pair is terminated at 
the respective local characteristic impedance 104 and 
the respective sampling arm 113 is connected to the 
second test channel receiver 112. The measurement 
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ports 103 of the direct pairs not being measured are 
also terminated in local characteristic impedances and 
the respective sampling arms 113 are terminated in 
respective local sampling arm terminating loads 116. 
The signal generator 105 is again swept through the 
desired frequency range and for each frequency point in 
the range, the VNA 200 measures a reflection response, 
a transmission response and a stimulus of the locally 
terminated thru 601 and stores the data in the 
following arrays : 

A f nm_ t e rmt h r u , 

Bf nm_termthru f and 

Rlfnm termthru. 



[28] With specific reference to FIGURE 9 of the 

drawings, and with the thru 601 still connected between 
the measurement ports 103 of the first direct pair, the 
transfer switch 106 is re-configured so that the signal 
generator 105 is in the second signal path 108 and the 
first signal path 107 is terminated in the 
characteristic impedance 109. The switch network 150 
is also reconfigured so that the measurement port 103 
in the first direct pair that is capable of connection 
to the second test channel receiver 112, measurement 
port 103 3 in the illustrated example, is connected to 
the second signal path 108 and the respective sampling 
arm 113 3 is connected to the second test channel 
receiver 112. The measurement port 103 in the first 
direct pair that is capable of connection to the first 
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test channel receiver 111, measurement port 103 x in the 
illustrated example , is terminated in the respective 
local characteristic impedance 104i and the respective 
sampling arm 113i is connected to the first test 
channel receiver 111. The measurement ports 103 of the 
direct pairs not being measured are locally terminated 
in their characteristic impedances 104. Additionally, 
the sampling arms 113 of the unused measurement ports 
103 are terminated in their respective local sampling 
arm terminating loads 116. The signal generator 105 is 
swept through the desired frequency range and for each 
frequency point in the range, the VNA 200 measures a 
reflection response of the locally terminated thru 601, 
a transmission response of the locally terminated thru 
601, and the level of the stimulus signal as measured 
at the second reference channel 202. The measured 
values are stored in data arrays: 

A r nm_termthru, 
B r nm_terrothru, and 
R2 rnm ^ termt h ru . 

[29] The same process is followed to measure and record 
similar measurements of the thru 601 for all remaining 
direct pairs and indirect pairs of the measurement 
ports 103. It is preferred to measure the same 
frequency points for all calibration measurements along 
the frequency range so that each array has a 
corresponding measured value for each frequency point. 
It is acceptable, however, to interpolate the data to 
obtain a value for a specific frequency value as long 
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as the frequency value is within the boundary of the 
lowest measured frequency and the highest measured 
frequency in the desired frequency range and the 
interval between measured frequencies is small enough 
to fully characterize the DUT including any resonances 
thereof. If multiple line calibration standards are 
used to obtain a broader frequency range, the 
measurements taken by the VNA 200 are stored in 
appropriate array elements in a larger array having an 
element for each frequency along the frequency range of 
interest. Accordingly, multiple steps of connecting 
calibration standards and making measurements may be 
performed to completely populate a single data array. 

[30] With specific reference to Figure 10 of the 

drawings, there is shown a flow graph between any first 
port and any second port of the VNA 200. A multiport 
embodiment has a similar calibration flow graph for 
each of the direct pairs. The calibration flow graph 
represents the forward directivity (shown as "EDF" in 
the f 511 patent and Sxn in the present teachings) 
1901, forward source match (shown as "ESF" in the '511 
patent and Sx 22 in the present teachings) 1902, and 
forward reflection tracking error coefficients (shown 
as "ETF" in the 1 511 patent and Sxi 2 Sx 2 i in the present 
teachings) 1903 for the X error adapter 1910. The 
calibration flow graph also represents reverse 
directivity (shown as "EDR" in the ! 511 patent and Sy u 
in the present teachings) 1904, reverse source match 
(shown as "ESR" in the '511 patent and Sy 22 in the 
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present teachings) 1905, and reverse reflection 
tracking (shown as "ERR" in the '511 patent and Syi 2 Sy 2 i 
in the present teachings) error coefficients 1906 for 
the Y error adapter 1920. Not shown on the flow graph, 
but four remaining error coefficients that are also 
used in the ten term error correction process are 
forward transmission tracking (shown as "ETF" in the 
'511 patent and inm in the present teachings where port 
"m" sources a stimulus and port "n" is the measurement 
port), forward load match (shown as "ELF" in the '511 
patent and r m in the present teachings for port "m" ) , 
reverse transmission tracking (shown as "ETR" in the 
f 511 patent and i mn in the present teachings where port 
"n" sources a stimulus and port "m" is the measurement 
port) , and reverse load match (shown as "ELR" in the 
f 511 patent and r n in the present teachings for port 
"n") . An embodiment of a method according to the 
present teachings determines the error coefficients for 
each direct and indirect pair based upon measurements 
of the TRL calibration standards for each direct and 
indirect measurement port pair. The flow graph 
represents an S-parameter matrix for the X error 
adapter 1910, S X / which corresponds to a subset of the 
error artifacts for a first measurement port 103 in the 
direct or indirect pair, and an S-parameter matrix for 
the Y error adapter 1920, S y , which corresponds to a 
subset of the error artifacts for a second measurement 
port 103 in the direct or indirect pair. The remaining 
error coefficients are calculated from the X and Y 
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error adapter error artifacts. 

[31] S-parameter matrix S act represents the S-parameters 
of a calibration standard without the contribution of 
the X and Y error adapters. The S-parameter matrix of 
the X error adapter may be expressed as T-parameters 
using the following known conversion where port 1 
represents a port from which a measurement is taken and 
port 2 represents a port that sources a stimulus 
signal : 



Tx Tx 
Tx 2i Tx 22 



1 


-Sx n 


SX 2\ 


Sx 2l 


Sx n 




Sx 2i 


Sx 2t 



(1) 



[32] Accordingly, the matrix Sx may be converted into 
corresponding T-parameters expressed as Tx. If the 
matrix T act _ thru expresses the T-parameters of just the 
thru 601 and T meas _thru expresses the T-parameters of the 
thru 601 as measured in context with the X and Y error 
adapters, then the following relationship holds true: 



T X T 



x 1 act_thruTy— T meas _thru 



(2) 



[33] Similarly, if the matrix T act _ii ne expresses the T- 
parameters of just the line 4 01 and T meas _iine expresses 
the T-parameters of the line 401 as measured in context 
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with the X and Y error adapters, then the following 
relationship holds true: 

T x T ac t_lineTy = T m eas_line ( 3 ) 

If the following relationships are defined: 

Tact_x ==r f'act_linel , act_thru ( 4 ) 

and 

Tmeas_x~T meas ^ii ne T meas ^thru ( 5 ) 

then the following equation can be written: 

TxT ac t_x == T n ieas_xTx ( 6 ) 

[34] The thru 601 and the line 401 are each assumed to 
be perfectly matched. Therefore, the value of their 
reflection coefficient in the respective actual S- 
parameter matrix is set to zero. If the thru 601 has a 
non-zero length transmission coefficient, it is defined 
by S2i_thru=Si2_thru - The line 4 01 has a transmission 
coefficient defined by S2i_iin e s= Si 2 _iine. From equation 
(4), therefore T act x may be expressed as: 
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s 2l 



T 

act x 



0 



0 

$2l thru 



(7) 



[35] Measurements of the source thru 601 and the source 
terminated line 401, each provide eight frequency 
dependent arrays of measured and stored results . 
There are four thru forward reflection and transmission 
arrays and four thru reverse reflection and 
transmission arrays. The arrays of measured data for 
the source terminated thru 601 are used in an 
algorithmic formulation in the S-parameter domain to 
compensate for the presence of the signal transfer 
switch 106 prior to calculation of the T mea s_x matrix. 
The effect of the transfer switch 106 of the 
uncorrected measured S-parameter for the thru and line 
calibration standards, S mea s_iine and Sm e as_thru, are 
corrected by the formulation given by: 
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J trans ferswttxh _ corrected 



( A, 



A. R2, 



R\ f R2 r R\, 



1- 



R2 .r RK 

R\, R2 r 



R\, R2 r R\ t 



f B, 



R2 .r RK 

Rl, R2 r 





A, Rl r 


R2, 


R\, R2 r 


1 - 


R2, R] r 




R\ f R2 r 


B r 


Bj RK 


R2, 


R], R2 r 



l- 



(8) 



where A / f B / ,R] / and /?2,are the forward direction raw 

measurement data f i.e. when the signal transfer switch 
106 is directing the signal generator 105 to the first 
signal path 107, and A r , B r , R\ r and R2 r are the reverse 
direction raw measurement data, i.e. when the signal 
transfer switch 106 is directing the signal generator 
105 to second signal path 108. After using this 
correction, the relevant S-parameter matrices are 

Smeasnm_thru_transferswitch_ corrected and 

Smeasi™_iine_transferswitch_corrected where "n" and "m" represent 
the port numbers between which the respective 
calibration standards are measured. _ 

[36] Using the first direct pair as an example, 
measurement ports 103i and 1033, & corrected S- 
parameter matrix of the thru 601 measured in cascaded 
combination with the X and Y error adapters for the 
first direct pair is expressed herein as 
S rae asi3_thru_transferswitch_corrected- The correction f ormulation 
shown in equation (8) uses the arrays: Afi3_ th ru, Bfi 3 _ th ru/ 
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Rlfl3_thru/ R2fi3_thru/ A r i3_ thru/ B r i3_ t hru/ Rlrl3_thru/ and 
R2 r i3_ t hru tO Calculate S me asl3_thru_transferswitch_corrected • 

Converting the S me asi3_thru_tranferswitch_corrected matrix to the 
corresponding T -parameters using equation (1) yields 

matrix T me asl3_thru_transferswitch_corrected • To obtain the 
S me a S 13_line_transferswitch_corrected matrix for the first direct 

pair, the correction formulation shown in equation (8) 
uses the arrays: Afi3_i ine / Bfi3_ii ne , Rlfi3_iine/ R2fi3_ii ne / 
A r i3_iine, B rl3 _i ine / Rl r i3_iine/ and R2 r i 3 _iine. Converting the 

Corrected S me asl3_line_transferswitch corrected matrix to the 

corresponding T-parameters, yields matrix 

T me asl3_line_transferswitch_corrected • The 

Tmeasl3_thru_transf erswitch_corrected and T meas i3_ii ne _transf erswitch_corrected 

matrices are used in equations (4) and (5) to calculate 
T a ct_x and T meas _ x . 

[37] Referring now to the general case, T x is the T- 

parameter matrix for the X error adapter and is defined 
by its matrix elements as: 



Tx = 



Tx u Tx n 

'1'X-^ | TXf^ 



(9) 



Tmeas x is also defined by its matrix elements, and is 
represented as: 



T T 

* meas _x\\ * meas jr 1 2 



T T 

x meas _x2\ meas _ x 



22 



(10) 
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From equation (5), T meas _ x for measurement ports 103i and 
103 3 , which is expressed as T mea si3_x, is calculated using 

the T me asl3_t:hru_transf erswitch_corrected and 

T me asi3_iine_transferswitch_corrected matrices . Accordingly: 
T - T T 

* meas 13 x meas 1 3 _ tine transfer-switch _ corrected * nteas 1 3 thru _ transferswiuh _ corrected 

[38] Using the relationship in equations (4), 

substituting the terms in equation (6), and eliminating 
the S 2 i_thru/S 2 i__iine term, the following general equation 
can be written: 



Tx 2 1 ( ~ T meas xU + ^4 T meas r , 2 T meas y 2 1 *+" ( T mtas r } i T meas x22 ) + T mea< . , r 22 ) (11) 



and 



7jC 22 ( ~T m(eas _x\\ + yj^meas _xl2 ? meas _x2\ + ( ^ meas _x\\ ^ meas _ x 22 ) ^ meas „ 1 22 ) (12) 
-* ^12 mtrtn _xl2 



Based upon the T-parameter to S-parameter conversion, 
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Tx Tx 

±±2i_ anc i — 22. in terms of the corresponding S-parameter 
Tx u Tx n 

error adapter matrix may also be expressed as: 



and 



T?2L = Sx u =B (13) 
Tx u " 



lZ»=Sx u - SXnSX2X =A (14) 
Tx n " Sx 22 



[39] As one of ordinary skill in the art can 

appreciate, equations (11) and (12) are equal. Because 
there is a square root in the solution, there are two 
possible mathematical solutions. The smaller valued 
solution, defined by B, corresponds to the directivity 
error coefficient, S x n 1901, of error adapter X. The 
larger valued solution, defined by A, is a mathematical 
combination of source match 1902 and reflection 
tracking 1903. 

[40] As mentioned before, at frequencies around 500 MHz 
and below, the dimensions of the line 401 become very 
large and not practical. Calculation of the directivity 
1901 and the solution represented by A for the lower 
frequencies uses the measurements taken of the two 
high-quality matched loads 501 instead of the line 601. 



31 



Agilent PDNO: 10021239 
Inventor: Vahe Adamian 



It is assumed that the matched loads 501 are perfectly 
matched to the measurement port and have a zero 
reflection coefficient. The same algorithmic 
formulations shown in equations (5) through (14) are 
used. To understand the usage of the measurement 
results from the two matched loads, note that the thru 
601 has a non-zero-length transmission coefficient 
defined by Si 2th ru=Si 2 thru • The matched loads 501 have a 
transmission isolation coefficient defined by 
S2iioad ==s i2ioad- Due to the high isolation between the 
matched loads 501, S 2 noad is close to a zero value. 
Accordingly, S 2 iioad is set to a very small, non-zero 
value, such as 10" 10 in order to avoid division by zero 
ambiguity in the S-parameter to T-parameter conversion. 
From this, T act _ x at the lower frequencies can be 
calculated and is given by: 



T 



io-'° 

0 



0 

io-'° 

^21 thru 



(15) 



[41] As shown before, by using equations (4) and (5), 
substituting results into equation (6) and eliminating 
the S 2 i__ t hru/10" 10 term, equations (11) and (12) are 
derived. The S-parameters from the matched loads 501 
are corrected using equation (8) to yield 

S me asl3_load_transferswitch_corrected/ which is then Converted 
USing equation (1) tO yield T me asl3_load_transferswitch_corrected - 
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The T me asl3_load_transferswitch_corrected term is USed tO Calculate 

T mea si3 x in place of the terms measuring the line 401. 
The calculations in equations (11) and (12), therefore, 
are the same as for the line 401. 
[42] A similar process is performed to calculate terms 
in the error adapter Y. Beginning with equations (2) 
and (3) and defining the following relationships: 

T aC t_y" ss Tact_thru" 1 T.ct.lii» ( 16 1 

and 

T m eas_y = Tmeas_thru T meas _ii n e (1*7) 

then the following equation can be written: 

Tact_yTy = T y T me as_y ( 18 , 

[43] With specific reference to Figure 10 of the 

drawings, the known conversion for T-parameter matrix 
for the error adapter Y in terms of the S-parameters 
where port 1 is on the right and port 2 is on the left 
when looking at the DUT 101 is: 
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Ty = 



\ 


-Sy„ 


Sy n 


Sy t2 


Sy 22 


Sy ]2 Sy 2i -Sy u Sy 2 2 


Sy l2 


Sy X2 



(19) 



[44] Accordingly, the matrix Sy may be converted into 
corresponding T-parameters expressed as Ty. The 

T me asl3_thru_corrected and T me asl3_line_corrected matrices have 

already been calculated and are used in equation (17) 
to calculate T mea si3_y/ where: 



1 meas 1 3 y 



T T 

1 meas 1 3 _ y II meas 1 3 _ y 1 2 

T T 

A meas 1 3 _ v 2 i A meas 1 3 _ y 2 2 



(20) 



[45] Using equation (18), substituting the relationships 
in equations (16) and (17) and eliminating the, 
S2ithru/S 2 iiine term, the following equation for the first 
direct pair can be written: 



Ty n ( "Tineas 13^1 1 + j/j La,» 1 3 _ y 1 2 T mhos \ 3 _ y 1 \ + ( T meas 1 3 _ y 1 1 ^mrav 1 3 _ v 22 ) Lfl» 1 3 _ y 22 ) j21) 

7y n 27^^,3^21 



and 
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Ty n (—Tmeas\3_y\ \ + ^^meas\3_ } \2^mtas \3_y2 1 + (^measlj _y\ I ^n K a.s\3_y22 ) + 1 3 _ ,i/22 ) ^2) 

Ty 2] 2r wi?avl3> . 2 , 



[46] From equation (19), ^-and in terms of the 

corresponding the S-parameters for the error adapter Y 
is also given by: 



and 



^H = -SV H =D (23) 



JhL=%A-- S y u =C (24) 
Ty lx Sy 22 



[47] As one of ordinary skill in the art can 

appreciate, equations (21) and (22) are equal and 
because of the square root have two solutions. The 
smaller value or the first solution, defined by Syn, 
corresponds to the directivity error of the Y error 
adapter. The larger value or the second solution, 
defined by C, corresponds to the error coefficient, 

Sy v Sy ^ -Sy u , for the Y error adapter. 

Sy 22 

[48] With reference to the portion of the calibration 
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procedure that measures the high reflect calibration 
standard as illustrated in FIGURE 2 of the drawings, 
the high reflect standard 301, is connected to one 
measurement port 103! of the first direct pair and the 
same high reflect standard 301 is disconnected from the 
measurement port 103i and is then connected to the 
other measurement port 103 3 of the first direct pair. 
With specific reference to FIGURE 10 of the drawings 
the following equation can be written: 



_ c Sx n Sx 2X T ac( _ reJha _ x (?5) 

^ meas _reflect x ~ «l I + , « r V ' 

I JA 22 1 at: i reflect _x 



where r meas _refiect_x is the measured reflection 

coefficient of the high reflect standard 301 at the 

measurement port capable of connection to the first 

test channel receiver, measurement port 103i in the 

first direct pair, and r ac t_refiect_x is the actual 

reflection coefficient of the high reflect standard at 

the same measurement port 103i. The same high reflect 

calibration standard 301 is connected to the opposite 

i . 
port in the first direct pair, measurement port 103 3 in 

the specific example. With respect to the error 

adapter Y, the following equation can also be written: 



Sy ]2 Sy 2i r actreJiecty (26) 

metis reflect _y s I I i rt r-» 

I *~>y 2 2 * act _ reflect _y 
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where r mea s_refiect_y is the measured reflection 
coefficient of the high reflect standard 301 at 
measurement port 103 3 and r act _ re f iect_y is the actual 
reflection coefficient of the high reflect standard at 
the measurement port 103 3 . A value for the measured 
reflection coefficient of the high reflect standard, 
r me as_refiect_x/ for the measurement port 103 x may be 
obtained from the following measured and stored arrays 
A re fiect_i/Rlrefiect_i. Similarly, a value for the measured 
reflection coefficient of the high reflect standard, 
r me as_refiect_ y , for the measurement port 103 3 may be 
obtained from the following measured and stored arrays 
Brefiect_3/R2 ref iect_3. Because the same high reflect 
standard is connected to measurement ports 103i and 
103 3 , it is possible to solve for r ac t_ref iect_x in 
equation (25) and r aC t_ref iect_ y in equation (26), and set 
the r terms -equal to each other. From the resulting 
relationship and equations (13), (14), (23), (24), (25) 
and (26), the following relationship can be written: 



meas _ reflect _ x 

)(C- 

^" meas _ reflect 

Sx 22 = . — - — ■ — _~ — = — l*'> 



( ^ ^nteas _ refket i _ x X ^ ^ meas _ reflect _ y ) 



As one of ordinary skill in the art appreciates, 
equation (27) has two unknown terms, but permits the 
expression of Sx 22 in terms of Sy 2 2 - Accordingly, 
another relationship is necessary in order to solve for 
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these two unknown terms . 

[49] With specific reference to FIGURES 7 and 10, the 
following equation can also be written: 



r _ s Sx n Sx 2] Sy 22 2Q 

1 ~ " \-Sx 22 Sy 22 



where r raea s_thruii is measured as A f 13 _ thru /Rlf i3_thru for the 
measurement port of the first direct pair that is 
capable of connection to the first test channel 
receiver 111. From (13), (14), (23) and (24), the 
following can be written: 



g | V~ x meas _ reflect _ x X ^ ^ meas _ reflect _ >» X ^ ^meas _ thru \ 1 ) (29) 



_ EES 



meas reflect _x 



X D ^mea.% _ re//ec / _ y X ^ ^meav _ /fcru I J ) 



Sx 2 2 is the source match error coefficient at the first 
measurement port 103. Due to the square root function 
in equation (29), there are 2 possible solutions for 
Sx 2 2- Additional information is needed, therefore, to 
solve for the appropriate solution for the present 
measurement . 

[50] The reflection-tracking error coefficient S x i 2 S x2 i 
is calculated for each of the two possible solutions to 
yield : 
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and 

S x 1 2 S x 21 — (B ~ A) S x22 solution 2 (31) 

[51] Squaring either equation (30) or (31) provides the 
complex value (S x i 2 S x2 i) 2 . Because the square of the 
reflection-tracking coefficient is a complex variable, 
its square root has two solutions. Each of the two 
solutions has the same magnitude, but they are 180 
degrees out of phase. Only one of the solutions is 
proper. To select the proper solution, both solutions 
for reflection tracking are solved at DC. Based upon 
the properties of reflection tracking, the argument of 
the proper solution is zero at DC. Given this boundary 
condition, therefore, the solution for reflection 
tracking with an argument closest to zero at a DC 
frequency is the proper solution and a value for the 
reflection tracking parameter may be determined. 

[52] Because S-parameter measurements are not taken at 
DC, the function that is fit to the argument of the 
reflection-tracking coefficient is extrapolated to a 
frequency value at DC. Specifically, the argument, 8, 
for each solution is calculated as a function of 
frequency. A first solution for the argument is 
represented as 8 ( f ) S x22_soiutioni and a second solution is 
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represented as 0 ( f ) S x22_soiution2 - For each solution, as 
the calculated function of frequency crosses through 
+180 or -180 degree points, 360 degrees is added to 
form a continuous function of phase as a function of 
frequency. Each resulting continuous function is fitted 
using a least sum squared fitting algorithm with a two 
term polynomial where f represents a frequency value 
and where: 



i 

0(f)s.x22 ^sotutionl = ^,4/' = 4) _sx 22 _.solution I + 4_sx 22 solution I / ( 32 ) 



r=0 

and 

i 



0( fXsx 22 _ solution 2 ~ ^ 4 f ~ 4_sx!!_ solution! + 4 _xx 22 _ solution! f (33) 



i=0 



[53] Each two term polynomial defines a straight line 
with slope, Ai_ sx22 , and y-intercept, A 0 _ S x22 • The y- 
intercept represents a value of 9 at DC. Of the two 
possible solutions for A 0 _ S x22/ only one of the solutions 
has a y-intercept of zero. Accordingly, whichever 
value, A 0 _ S x22_soiutioni or A 0 _ S x22_soiution2f is closest to a 
zero value is the proper solution. When the proper 
solution is identified, both the magnitude and phase 
are then calculated for the reflection tracking 
coefficient, S xl2 S x2 i. From the calculation of the 
reflection tracking coefficient, source match error, 
S x2 2/ is calculated where: 
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c _ Sxi2§x2\ ( 34 ) 

' V22 ~ B-A 

The factors represented as "B" and "A" are known from 
equations (13) and (14) . Accordingly, the source match 
error, S x22 / is also fully determinable. 
[54] When a value for Sx 22 is known, a value for Sy 2 2 may 
be calculated from equation (27) . Sy 22 is the source 
match error coefficient of error adapter Y at the 
measurement port capable of connection to the second 
test channel receiver 112, which is measurement port 
103 3 in the specific example. 

[55] A definite value for Sy 22 and equations (23), (24) 
and (27) permit calculation of a reflection tracking 
for the error adapter Y and is given by: 

Sy l2 Sy 2i =(D-C)Sy 22 ( 35 ) 

[56] Accordingly, the directivity, source match and 
reflection tracking for the error adapter X and the 
error adapter Y for the first direct pair is 
determined. The same process is followed to calculate 
the directivity, source match, and reflection tracking 
for the X and Y error adapters for each direct and 
indirect pair. The X error adapters are defined as the 
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error artifacts presented in series with the 
measurement ports 103 capable of connection to the 
first test channel receiver 111. Similarly, the Y 
error adapters are defined as the error artifacts 
presented in series with the measurement ports 103 
capable of connection to the second test channel 
receiver 112 . 

[57] In the specific four-port embodiment, the 

measurements and calculations described herein for 
measurement ports 103i and 103 3 yield directivity, 
source match and reflection tracking for the error 
adapter X related to measurement port 103i and 
directivity, source match, and reflection tracking for 
the error adapter Y related to measurement port 103 3 . 
The same measurements and calculations described herein 
for measurement ports 103 x and 103 3 are performed for 
the second direct pair and if source terminated 
measurements are taken for the indirect pairs, all 
indirect pairs as well. If source terminated 
measurements are not taken for the indirect pairs, 
directivity, source match, and reflection tracking are 
determined according to the teachings US patent 
application Ser. No. 09/954,962 (herein "the '962 
patent application"), the teachings of which are hereby 
incorporated by reference. In a certain multiport 
embodiment, the same measurements and calculations are 
made for each direct pair and indirect pair to yield 
directivity, source match and reflection tracking for 
the error adapter X related to the measurement port of 
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the direct or indirect pair capable of connection to 
the first test channel receiver 111 and directivity, 
source match, and reflection tracking for the error 
adapter Y related to the measurement port of the direct 
or indirect pair capable of connection to the second 
test channel receiver 112. Accordingly, an N port DUT 
100 has N(N-l)/2 different X and Y error adapters 
associated therewith. In another multiport embodiment, 
the disclosed measurements are made for all direct 
pairs to calculate directivity, source match, and 
reflection tracking and the error terms for the 
indirect pairs are calculated according to the f 962 
Patent Application. 

[58] Using the forward reflection and transmission 
measurements made on the locally terminated thru 601, 
which in a specific example of the first direct pair 
are the A f i3_termtnru, Bfi3_termthru/ and Rlf i3_termthru arrays, it 
is possible to solve for a load match error coefficient 
presented at the measurement port capable of connection 
to the second test channel receiver 112 and a forward 
transmission tracking error coefficient for the first 
direct pair. The load match for the measurement port 
103 3 , r /3 , and forward transmission tracking for the 
first direct pair,r l3 , are given by: 



j_ _ $ X l\ icrmthrv ^ ^\ f }3 ^urmthru) (36) 

Sx u Sx 22 — Sx x2 Sx 2l ~Sx 22 (A fnierm(hrv /R] n3lvrntthru ) 
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and 



r„ =(B J , Klemthn /R\ /rj Kmlhn )(\ -& a r u ) ( 37 ) 

'59] Using the reverse reflection and transmission 
measurements made on the locally terminated thru 601, 
which in a specific example of the first direct pair 
are the A r i3_tennthru/ B ri 3_termthru/ and R2 f i3_termthru arrays, it 
is possible to solve for a load match error coefficient 
presented at the measurement port capable of connection 
to the first test channel receiver 111 and a reverse 
transmission tracking error coefficient for the first 
direct pair. The load match for the measurement port 
103i, r n/ and the reverse transmission tracking 
coefficient, r 3)/ are given by: 



j_ ^V)| ( Br I 3 _ termthru ^ ^ r 1 3 _ termthru ) (38) 

1,1 Sy u Sy 22 -Sy i2 Sy 2i -Sy 22 (B rl3 lemthru I R2 rl3 lermlhru ) 



and 



[60] In the specific illustrated embodiment, using the 
measurements made of the locally terminated thru 601 
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for the second direct pair, the same algorithmic 
formulations shown in equations (36) through (39) as 
described for the first direct pair are applied to the 
measurement ports of the second direct pair. 
Accordingly, forward and reverse directivity, source 
match, reflection tracking and load match error 
coefficients for each measurement port in the direct 
and indirect pairs and forward and reverse transmission 
tracking error coefficients for all direct and indirect 
pairs are determined . 

[61] In a multiport embodiment of a method according to 
the teachings of the present invention, the forward and 
reverse directivity, forward and reverse source match, 
forward and reverse reflection tracking, forward and 
reverse load match, and forward and reverse 
transmission tracking error coefficients for each 
measurement port in all of the direct pairs are 
similarly determined . 

[62] Because the directivity, source match and 
reflection tracking are fully determinable, 
characteristics of the high reflect calibration 
standard 301 may also be calculated. Using equations 
(25) or (2 6) and solving for r act reflect x or T a 

ct reflect yr 

yields : 



r — r 

r — L meas ^reflect _x { A f) \ 

1 act effect _x ~ /j_ r {HU} 
V* 1 mea.s _ reflect _ x / ° x22 
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or 



n — r 

'-' 1 nieas _ reflect y 



^act _ reflect _ V t^> _ p \ O 

W *■ meas rejlea _y ) y22 



[63] . The argument of r act _refiect_x as a function of 

frequency, where 360 degrees is added as the function 
passes through the +/-180 degree points as previously 
described, is fitted to a least sum squared fitting 
algorithm with a two term polynomial as in equations 
(32) or (33) . The resulting two term polynomial that 
describes the proper function includes Ao_ract_refiect_x as 
the y-intercept and Ai_ ra ct_ref iect_x as the slope. If 
A 0 ract reflects h'as a value closest to 180 degrees, then 
the high reflect calibration standard 301 has a short 
circuit reflection coefficient. If Ao_ract_refiect_x has a 
value closest to 0 degrees, then the high reflect 
calibration standard 301 has an open circuit reflection 
coefficient. Accordingly, the type of high reflect 
calibration standard measured can be deduced from the 
least sum square fit calculations. If the electrical 
length of the thru 601 is unknown, however, it is 
necessary to know a phase offset of the high reflect 
calibration standard 301 for additional correction. If 
t he electrical length of the thru is known, it is 
possible to determine a type of high reflect 301 and 
calculate the phase offset of the reflect 301. 

[64] The forward and reverse transmission tracking 



Agilent PDNO: 10021239 
Inventor: Vahe Adamian 



error coefficients of the proximal pairs are: r l2 , r 2 , , 
r 34 , and r 43 , which may be determined either through 
measurement and calculation or through pure 
calculation. The method for pure calculation of the 
forward and reverse transmission tracking error 
coefficients is taught in US Patent Application Ser. 
No. 10/098,040, the teachings of which are hereby 
incorporated by reference. 

[65] Alternatively and with specific reference to 
FIGURE 11 of the drawings, the forward transmission 
tracking error coefficient for each proximal pair where 
both measurement ports 103 of the proximal pair are 
capable of connection to the first test channel 
receiver 111 is measured and calculated by connecting 
the thru 601 between the measurement ports 103 of the 
proximal pair, connecting the signal generator 105 to a 
first one of the measurement ports 103 of the proximal 
pair and terminating the sampling arm 113 in a local 
terminating impedance 116. For purposes of 
nomenclature, the first measurement port 103 of the 
proximal pair is designated as "port F" . The other 
measurement port 103 of the proximal pair, which for 
purposes of nomenclature is designated as "port G", is 
terminated in a local terminating impedance 104, and 
the respective sampling arm 113 is connected to the 
first test channel receiver 111. The VNA 200 measures 
and stores a ratio of a transmission response over the 
reference signal, Af F G_termthru/RlfFG_termthru- The ratio is 
used in the transmission tracking error coefficient 
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equation where: 

r WV = i^JFG jermthru ^ ^ JhXi jermihru) 0 ~ $ X 22 _ port/ ^ portl-} (41) 

[66] The reverse transmission tracking error 

coefficient for the same proximal pair, ports F&G, 
where both measurement ports 103 of the proximal pair 
are capable of connection to the first test channel 
receiver 111 is measured and calculated by keeping the 
connection of the thru 601 between the measurement 
ports 103 of the proximal pair, connecting the signal 
generator 105 to the second one of the measurement 
ports 103 of the proximal pair and terminating the 
sampling arm 113 in a local terminating impedance 116. 
The first one of the measurement ports 103 of the 
proximal pair is terminated in the local terminating 
impedance 104, and the respective sampling arm 113 is 
connected to the first test channel receiver 111. The 
VNA 200 measures and stores a ratio of a transmission 
response over the reference signal, 
A rFG _termthru/RlrFG^termthru. The ratio is used in the 
transmission tracking error coefficient equation where: 

T - i^rhd Jermthrv ^ ^ rhXi jermthru) 1 0 " ^22 _ ponO ^porH i ) (42) 

[67] For purposes of illustration and with 
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specific reference to FIGURE 11, there is shown 
connection diagrams for the determination of the 
forward and reverse transmission tracking error 
coefficients for the proximal pair comprising 
measurement ports 103 x and 103 2 . The forward 
transmission measurement is made by connecting the 
signal generator 105 to the first signal path 107. The 
switch network 150 is configured so that the first 
signal path 107 is connected to the measurement port 
103iand the respective sampling arm 113i is terminated 
in the local sampling arm impedance 116i. The switch 
network 150 is further configured so that the 
measurement port 103 2 is terminated in the local 
terminating impedance 104 2 and the sampling arm 113 2 is 
connected to the first test channel receiver 111. The 
signal generator 105 sweeps through the plurality of 
frequencies that define the desired frequency range and 
measures the ratio A f i2_termthru/Rlfi2_termthru - Using 
equation (40) the forward transmission tracking error 
coefficient for the proximal pair is calculated as: 

2 = ( ^/ 1 2 _ itrmthru ^ ^ J 1 2 _ lermthrv ) ( ' ~~ ^*22 _ part I ^ purl 1 ) 

[68] The connection of the thru 601 between the 
proximal pair and the signal transfer switch 106 
configuration are maintained. With specific reference 
to FIGURE 12 of the drawings, the switch network 150 is 
reconfigured so that the measurement port 103 2 is 
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connected to the first signal path 107 and the 
respective sampling arm 113 2 is terminated in the local 
sampling terminating impedance 116 2 . Additionally/ the 
switch network 150 is configured so that the 
measurement port 103i is terminated in the local 
terminating impedance 104i and the sampling arm 113i is 
connected to the first test channel receiver 111. The 
signal generator 105 sweeps through the plurality of 
frequencies that define the desired frequency range and 
measures the ratio A r i 2 _termthru/Rlri2_termthru - Using 
equation (41) the reverse transmission tracking error 
coefficient for the proximal pair is calculated as: 

**2 1 = ( A r 1 2 _ mrntihru ^ ^ ' r 1 2 lermthrv )( ' ~~ ^ 22 _ pori 2 ^port 2 ) 

[69] The same measurement and calculation process is 
repeated for the remaining proximal pairs, which in the 
illustrated embodiment is the proximal pair comprising 
measurement ports 103 3 and 103 4 . As one of ordinary 
skill in the art appreciates, the measurements for the 
proximal pair capable of connection to the second test 
channel receiver 112 are carried out using the same 
process, but with the second test channel receiver 112 
and the second reference channel 202 as measurement 
devices. In a multiport embodiment, the measurement 
and calculation process is repeated for all of the 
proximal pairs . 

[70] When all ten error coefficients are known for each 
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direct, indirect, and proximal pair of measurement 
ports, the DUT is inserted and measured by the VNA. 
The resulting measured S-parameters for each direct and 
indirect pair of the inserted DUT are corrected 
according to the '511 patent to arrive at actual S- 
parameters for the DUT, Sactii/ S ac t22/ S a cti2/ and S ac t2i for 
each direct and indirect pair. U.S. Patent Application 
Ser. No. 09/954,962 (herein "the '962 Patent 
Application"), the teachings of which are hereby 
incorporated by reference, teaches a process for 
accepting a plurality of actual S-parameter matrices, 
each actual S-parameter matrix representing one direct 
or indirect pair. Multiple direct and indirect pairs 
are used to represent a single multiple-port DUT. 
According to the teachings of the '962 Patent 
Application, the multiple actual S-parameters for the 
two-port devices that combine to represent the multi- 
port DUT are normalized to impedances presented to the 
measured ports by the multiport test set. The multiple 
S-parameter matrices are combined into a larger 
normalized single S-parameter matrix to represent the 
corrected and normalized S-parameters that represent 
the DUT as a whole. The larger S-parameter matrix is 
further normalized to a characteristic impedance, for 
example 50ohms. For purposes of the remaining 
disclosure, the measured, corrected and normalized S- 
parameter matrix for the DUT is referred to as the 
matrix S. 

[71] Because the calibration process uses a non-zero 
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through calibration standard 601 and perhaps an offset 
high reflect calibration standard, a reference plane 
based upon the through calibration standard 601, termed 
a "shifted calibration reference plane", is placed in 
the center of the electrical length representing the 
aggregate of the electrical lengths of any delay of the 
non-zero through calibration standard 601 and any 
offset of the high reflect calibration standard 301. 
In on-wafer measurements, it is typical to have a zero 
offset high reflect calibration standard and an unknown 
non-zero through calibration standard. For coaxial 
measurements, it is possible to calibrate with a zero 
delay through and an offset high reflect. The 
teachings herein accommodate both situations as well as 
other combinations of through and high reflect 
calibration standards not specifically mentioned. 
FIGURE 6 of the drawings, illustrates a placement of 
the shifted calibration reference plane 701. The 
proper reference plane for measurement of the DUT, 
however, is a reference plane that is disposed at each 
measurement port that connects to the DUT, termed "the 
DUT reference plane". The DUT reference plane for port 
1 is shown as 7 02 in FIGURE 6 of the drawings. Each 
measurement port 103 has a respective and similarly 
situated DUT reference plane. * In order to establish 
the DUT reference plane for all measurements, S- 
parameters for each direct and indirect pair are phase 
shifted. The amount of phase shift required for the S- 
parameters is calculated from a determination of the 
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electrical length of the non-zero through calibration 
standard. The mispositioning of the reference plane may 
be corrected by determining the electrical length of 
the combination of the thru 601 and the offset high 
reflect calibration standard 301 and adjusting 
respective measured and corrected S-parameters 
accordingly . 

[72] If the electrical length of the thru 601 is 
unknown and the offset of the high reflect 301 is 
known, the electrical length of the non-zero thru 601 
between each direct and indirect pair, may be 
determined. Begin with the least sum squared fit 
(refer to equation (32) and accompanying text) of the 
measured reflection parameter r act _ re fiect_x as an example 
(refer to equation (40) and accompanying text) of the 
measurement port for the direct or indirect pair of 
interest. The slope of the function, Ai_ ra ct_ref lectx* 
represents a derivative of the argument of the high 
reflect calibration standard 301. The derivative of 
the argument defines a change in phase with respect to 
the corresponding change in frequency as a function of 

frequency or — . Any phase delay as a function of 

frequency that is measured when the high reflect 
calibration standard 301 is connected that cannot be 
attributed to the high reflect calibration standard 
offset can be attributed to the phase delay introduced 
into the measurement as a result of the calibration 
process using the non-zero through 601. If the 
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electrical length of the non-zero thru is known, then 
any phase delay that cannot be attributed to the thru 
calibration standard 601 can be attributed to the 
offset of the high reflect calibration standard 301. 
This phase delay attributable to the non-zero thru 601 
and/or the offset reflect calibration standards, 
therefore, may be removed from the final DUT 
measurement when the electrical lengths are known. In 
the case of a zero offset high reflect standard, the 
phase shift of the high reflect standard 301 as a 
function of frequency is equal to the negative of the 
phase shift introduced by the non-zero thru 601. The 
aggregate electrical delay of the non-zero through and 
the offset of the high reflect calibration standard, 
represented as T in seconds, is defined as: 



Accordingly, the electrical delay, T, as a function of 
frequency is expressed as: 



T = - 



(43) 



360^/ 



T = 




(44) 



360 



[73] The electrical length, 1, in centimeters is 

calculated from the electrical delay as a function of 
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frequency where c is equal to 2.99793xl0 10 or the speed 
of light in centimeters per second and: 

I = cT (45) 



The process to determine the electrical length, 1, of 
the aggregate offset from the non-zero thru 601 and/or 
offset high reflect is repeated for each one of the 
direct and indirect pairs. The electrical length may 
be calculated or measured for the indirect and proximal 
pairs. As a practical matter in the case of on-wafer 
measurements, even if the same on-wafer calibration 
standard is used, a different landing is made for each 
direct and indirect pair. Because the through 
calibration measurement for each direct and indirect 
pair has a different landing, the electrical length of 
the through calibration standard will be different for 
each direct and indirect pair. Accordingly, a 
different electrical length is maintained for each 
direct, indirect, and proximal pair. 

[74] In many cases, a thru 601 between proximal pairs 
is not commercially available on-wafer for use with 
differential probes. Accordingly, the electrical 
length of the on-wafer thru 601 for proximal pairs is 
calculated instead of measured for purely practical 
reasons. If an on-wafer thru calibration standard for 
use with differential probes were to become available, 
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the electrical length between proximal pairs would be 
measured and calculated as it is for the direct and 
indirect pairs of measurement ports. Until then, the 
electrical length is calculated from an average of the 
direct and indirect pair having the proximal pair 
measurement ports in common. 

[75] In an embodiment where a thru 601 may be measured 
between proximal pairs, for example coaxial 
measurements, the electrical length of the thru 601 
between proximal pairs may be measured and calculated 
as described for the direct and indirect pairs. 

[76] In an embodiment where it is not possible or 
practical to measures a thru 601 between proximal 
pairs, the electrical length of the through calibration 
standard 601 is measured for the direct and indirect 
pairs. The electrical length of each proximal pair is 
then calculated from an average of the direct pair and 
the indirect pair having only one port in common. As 
an example, the electrical length of the thru 601 
between ports 3 and 4 is the average of the electrical 
length of the thru calibration standards between ports 
1 and 3 and between ports 1 and 4 . 

[77] For purposes of nomenclature, electrical length is 
represented as "1" with a subscript that represents the 
two ports to which it applies. Accordingly, the 
electrical length between measurement ports 103i and 
103 3 is represented as "li 3 " • In the illustrated four 
port DUT, electrical length is determined for all 
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through calibration standards yielding l i3/ 1 2 4/ li4/ and 
I23 electrical lengths for the thrus 601 measured 
between the direct and indirect pairs. The electrical 
length calculation of the two remaining proximal pairs 
yields I12 and I34. 

[78] In order to make the appropriate calculations to 
yield the electrical length of the direct and indirect 
pairs, the measurements made for the source terminated 
thru are used. Specifically, for the thru 601 measured 
between measurement ports 103i and 103 3 , the 
measurement arrays Afi3_ t hru, Rlfi3_tnru, Bfi3_ t hru/ R2f i3 _ thru , 
A r i3__thru, Rlri3_thru, B r i 3 _ thru , and R2 r i 3 _ thru are used to 
calculate the r ac t_ref iect_x term for the first direct 
pair. The r ac t_ref iect_x is used to obtain the 
Air_act_refiect_x term from equation (32) or (33) to make 
the electrical length determination from equations 
(41)-(43). A similar calculation is made for all 
direct and indirect pairs of the multiport 
configuration to determine the electrical length of 
each combination of the through calibration standard 
601 and offset reflect. 

[79] The elements of the matrix S are phase shifted 
according to the measured electrical length of the 
specific thru 601 and high reflect 301 used in the 
calibration process to correct the mispositioning of 
the calibration reference plane. For example, in a 
four port DUT embodiment, the electrical length l i3 is 
used to correct parameters Sn, Si 3 , S 3i , and S 33 . 
Similarly the electrical length 1 24 is used to correct 
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parameters S 22 , S 24 , S 42 , and S 44 , electrical length l i4 is 
used to correct parameters Si 4/ and S 4i , electrical 
length 1 2 3 is used to correct parameters S23/ and S32/ 
electrical length 1* 2 is used to correct parameters Si 2/ 
and S 2 i, and electrical length 1 34 is used to correct 
parameters S 34 , and S43. It is assumed that the same 
high reflect calibration standard 301 is used for all 
measurement ports 103. 

[801 In a currently preferred method, the resulting 
measured, corrected, and normalized S-parameters are 
adjusted at each frequency. If the unrotated S- 
parameter is represented as a general complex 
parameter, then: 

S =\p\e~'*> (46) 

where 9 0 is the phase of the corrected and unrotated S- 
parameter. The 9 term may be represented as a 
component of measured phase attributable to the 
measured DUT, 9 du t/ and a component of measured phase 
attributable to the shifted reference plane, 69, where: 

o**=o* + se (47) 

The length, 1, in centimeters is then used to calculate 
phase shift attributable to the non-zero thru 
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calibration standard 601 and high reflect standard as a 
function frequency, f, in GHz where: 

S0(f) = -12.0083// (48) 

where 69 is determined as a function of frequency for 
each electrical length represents the phase rotation 
attributable to the reference plane shift. When the 
electrical length and, therefore, the phase shift 
attributable to the sum of the electrical delay of the 
non-zero through calibration standard plus any offset 
in the high reflect standard as a function of frequency 
for each two-port combination is known from equation 
(49), the corrected and unrotated S-parameters may be 
adjusted as a function of frequency of the reference 
plane for each S-parameter . 
[81] Accordingly, phase corrected S-parameters are 
calculated as a function of frequency where: 

The resulting s-parameters in the S dut matrix represents 
measured, corrected and normalized S-parameters for the 
measured DUT where the measurement reference plane has 
been correctly placed at the DUT measurement reference 
plane 702 . 
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[82] In an alternative embodiment of an electrical 
length calculation for a TRL calibration process, 
forward and reverse measurements of the high reflect 
301, line 501 and optional matched loads, and source 
terminated through 601 calibration standards are taken 
for just the direct pairs. Directivity, reflection 
tracking and source match as well as electrical length 
of the thru 601 plus any offset high reflect 301 is 
measured and calculated for the direct pairs of 
measurement ports 103. Forward and reverse 
measurements of the locally terminated thru calibration 
standard 601 are then taken for just the direct and 
indirect pairs of measurement ports and calculated as 
taught herein. Alternatively, only the direct pair is 
measured and transmission tracking and load match may 
be calculated for the indirect and proximal pairs from 
the teachings of the '04 0 Patent Application. 
Optionally, forward and reverse measurements of the 
locally terminated thru 601 may be taken for 
calculation of transmission tracking and load match 
according to the teachings herein. The S-parameters of 
the DUT are measured, corrected and normalized 
according to the teachings herein. From the error 
terms, the electrical length of the thru plus any 
offset high reflect may be calculated for the indirect 
and proximal pairs as follows where "n" is a first port 
and "m" is a second port. Begin with the 
relationships : 
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r =EDF + &Ajonn (50) 

* SM portn portn . /TO IT p 

I - aor^,, I sA_ponn 



where T S m portn is the measured source match for port 
"n", r SA _ P ortn is the actual source match for port "n" , 
and EDFportn, ERF portn/ ESF portn are the calculated error 
coefficients for port "n" and: 



e- c r 

^ ° 2 1 _ thru _ nm 12 //trw ww LA_portm (51) 

SA _ purity H_ f/iru _ " w ~^ 1 O p 

1 ^ 22 _ thru _nm 1 LA_portm 



where r LA _ p0 rtm is the actual load match for port "m" and 
Sn_thru_nm/ S 2 i_ t hru_nmS i 2 _thru_run and S 22 _thru^nm are the measured 
S-parameters of the source terminated thru calibration 
standard 601. Because the thru calibration standard 
601 is defined as having a perfect match, equation (51) 
reduces to: 



c c p 

_ n Lj 2]_fhru nm Lj \2_ihru_nm _ LA_portin _ O O p (^9) 

1 SA por( „ = U + - j ^2\_thrv wm°12 thru nm 1 LA_portm v u * 1 



and 



L SA_ portn _ o o ( 53 ) 

~Z, — °2\ _ihru _ mv Lj \2_lhru_nm X ' 

* LA _ portm 
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Because S 2 i_i;hru_nm = Si2_thru_nm/ there are two possible 
solutions. Fitting the argument of the square root of 
S 2 i_thru_nm to a straight line using a least sum squared 
fit and equation (32) or (33), the solution with a y- 
intercept closest to zero is the correct solution and 
the resulting argument of the correct solution is the 
length of the electrical delay between the two proximal 
pair measurement ports. This length is used to correct 
for the shifted reference plane according to equations 
(46) -(49) and the accompanying text. 

[83] With specific reference to FIGURES 13 through 17 
of the drawings, there is shown a flow chart of an 
embodiment of a method according to the present 
teachings in which measurements for the direct and 
indirect pairs of measurement ports 103 are taken for 
the high reflect calibration standard 301, the line 
calibration standard 401, the source terminated thru 
calibration standard 601 and the locally terminated 
thru calibration standard 601. A measurement of the 
locally terminated thru calibration standard 601 is 
taken for each proximal pairs of the measurement ports 
103. The high reflect calibration standard 301 is 
connected 2201 to a measurement port 103 and the switch 
network 150 is configured for measurement 2202 by the 
VNA 200 of a ratio of the reflection response over the 
stimulus. See FIGURE 2 of the drawings. The ratio 
yields a value for a reflection ratio for a number of 
frequencies in a desired frequency range. The numbers 
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are stored in a data array where each element of the 
data array holds the measured ratio at a single 
frequency. The reflect calibration standard 301 is 
then disconnected and reconnected 2203 to a next 
measurement port, the switch network 150 is 
reconfigured, and similar measurements are made and 
stored in respective data arrays. The process 
continues for each measurement port 103 of the VNA 200. 
In a preferred embodiment, the desired frequency range 
for which all measurements are taken is the same. In 
this case, each element in the data arrays represents 
measured results at the same frequency point along the 
desired frequency range. 
[84] The flow chart of FIGURE 13 continues with the 

step of connecting 2301 and configuring 2302 the switch 
network 150 for termination of the line 401 at the 
source 109 and measurement of the line calibration 
standard 401 between a first direct or indirect pair. 
The VNA 200 measures 2302 a forward direction 
reflection and transmission response at the measurement 
ports 103 of the direct or indirect pair at the first 
and second test channel receivers 111, 112 as well as 
the first and second reference channels 201, 202. The 
switch network 150 is then reconfigured 2303 for the 
reverse direction measurement and the VNA 2 00 then 
measures 2303 a reverse direction reflection and 
transmission response at the measurement ports 103 of 
the direct pair at the first and second test channel 
receivers 111, 112 as well as the first and second 
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reference channels 201, 202. Not shown in the flow 
chart is the connection and measurement of the matched 
loads 501 as shown in FIGURE 6 of the drawings to 
extend the calibration to the lower frequency range. A 
measurement of the line 401 calibration standard is 
performed for each direct and indirect pair. 

[85] With specific reference to FIGURES 14 of the 
drawings, the flow chart continues with the step of 
connecting 2401 the thru 601 to a first direct pair of 
measurement ports 103. The switch network 150 is 
configured 2402 for forward direction measurement of 
the thru 601 with a termination at the source 109, and 
the forward direction reflection and transmission 
responses and the reference channel signals are 
measured and stored in data arrays. The switch network 
150 is then reconfigured 2403 for the reverse direction 
measurements with the thru 601 terminated at the source 
109 and the reverse direction reflection and 
transmission responses and the reference channel 
signals are measured and stored in data arrays. 

[86] The thru 601 remains connected and the switch 
network 150 is reconfigured 2501 for a forward 
direction measurement, where the thru 601 is locally 
terminated in a local impedance 104 within the switch 
network 150. The forward direction reflection and 
transmission responses of the locally terminated thru 
601 are measured and stored as well as the reference 
channel signals. The switch network 150 is then 
reconfigured 2502 for a reverse direction measurement 
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of the locally terminated thru 601, the reverse 
reflection and transmission responses and the reference 
channel signals are measured and stored. The process 
repeats 2503 for all N/2 direct pairs and N(N-2)/4 
indirect pairs of measurement ports 103 for an N port 
DUT . 

[87] At this point in the process, there is sufficient 
information to calculate an electrical length of the 
aggregate thru/offset reflect between direct and 
indirect pairs of measurement ports 103 used in the 
calibration process. For purposes of process, however, 
all calculations are made after all measurements are 
taken . 

[88] With specific reference to FIGURE 15, there is 
shown an optional portion of the measurement process 
wherein forward and reverse measurements of the locally 
terminated thru calibration standard 601 is taken 
between each proximal pair of measurement ports 103. A 
user connects the thru 601 between one of the proximal 
pairs 2504. The switch network 150 is then configured 
2505 for a forward measurement of the locally 
terminated through, measurements are taken and stored. 
The switch network 150 is then configured 2506 for the 
reverse measurement of the locally terminated thru 601, 
measurements are taken and stored. The process repeats 
for each N(N-2)/4 proximal pairs. 

[89] With specific reference to FIGURE 16 of the 

drawings, the flow chart continues with the steps of 
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calculating 2705 a forward reflection tracking error 
coefficient 1903 based upon a boundary condition at DC 
and using the forward reflection tracking error 
coefficient result, calculating 2706 a forward source 
match 1902 and forward directivity 1901 error 
coefficients. The reverse reflection tracking error 
coefficient 1920 is calculated 2707 also based upon a 
boundary condition at DC. Using the reverse reflection 
tracking error coefficient 1906 result, the reverse 
directivity 1904 and reverse source match 1905 are 
determined 2708. The method then determines 2709 the 
forward and reverse load match and the forward and 
reverse transmission tracking error coefficient to 
arrive at ten error terms based upon the TRL 
calibration. All ten error terms are determined for 
all direct and indirect pairs in a loop 2710 that 
increments 2711 through the plurality of two-port 
measurements . 
[90] With specific reference to FIGURE 17 of the 

drawings, the process continues with optional process 
paths shown as alternative paths on the flow chart. In 
one embodiment, the forward and reverse transmission 
tracking error coefficients for each proximal pair are 
calculated 2803 according to the teachings of the '040 
Patent Application. In an alternate embodiment, if the 
locally terminated thru 601 is measured 2804 between 
each proximal pair of measurement ports 103. The 
forward and reverse transmission tracking coefficient 
of each proximal pair are determined 2805 from the 
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measurements made. When the transmission tracking 
error coefficients are determined for each proximal 
pair of measurement ports 103, the process continues 
with the steps of determining 2801 a type of high 
reflect calibration standard 301 used in the 
calibration process and calculating 2802 an electrical 
length of each thru calibration standard 601 used in 
the calibration process between each direct and 
indirect pair. An electrical length is then typically 
calculated 2806 for each proximal pair based upon an 
average of the electrical length between direct and 
indirect pairs having a single common measurement port. 
If and when it is practical to measure a source 
terminated thru 601 between proximal pairs, the 
electrical length may be calculated as for the direct 
pairs . 

[91] With specific reference to FIGURE 18 of the 

drawings, when all of the systematic error coefficients 
are determined, the DUT 101 is inserted for measurement 
2902 to yield N(N-l)/2 2-port S-parameter matrices to 
represent an N port DUT. The measured DUT data is then 
corrected 2903 according to the teachings of the '511 
Patent. As also taught in the 1 511 patent, the 
corrected S-parameters matrices representing the 
plurality of 2-port devices are combined to yield an 
NxN S-parameter matrix that is corrected for the 
systematic errors of the VNA. The corrected S- 
parameter matrix is then normalized to a characteristic 
impedance . 
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[92] Using the respective electrical length determined 
for each thru 601 the NxN corrected and normalized S- 
parameters of the measured DUT are further corrected by 
shifting the reference plane from the shifted 
calibration reference plane 601 to the DUT reference 
plane 702 as a function of frequency. The reference 
plane shift operation yields a NxN Sf ina i matrix, which 
represents the S-parameters for only the DUT after all 
corrections are made. 

[93] In yet another alternative embodiment of a process 
according to the present teachings, reflection 
measurements for all measurement ports are taken, and 
measurements are taken for the line and the source 
terminated through for the direct pairs only. At this 
point in the process, it is possible to calculate 
electrical length for the thru/offset between the 
direct pairs. Additional measurements of the locally 
terminated thru 601 are taken for the direct and 
indirect pairs of measurement ports 103. The 
transmission tracking and load match coefficients may 
be calculated from the terminated thru measurements. 
The characteristics of the high reflect standard are 
determined and the electrical length for the indirect 
pairs is calculated based upon equations (50) - (53) . 
The electrical length for the proximal pairs may then 
be calculated based upon the average of the direct and 
indirect pair having a common port. All error 
coefficients and electrical lengths being determined, 
the DUT is measured, corrected, normalized and a 
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reference plane is shifted to a DUT reference plane as 
appropriate based upon the respective electrical 
lengths . 

[94] A method according to the present teachings may be 
performed by a computer with access to measured S- 
parameter characterization data for the through- 
reflect-line calibration standards as well as the S- 
parameter characterization data for the DUT. The VNA 
performs the measurements and either stores the 
resulting data arrays on removable media or may be 
connected to a computer via a communications bus for 
storage on some other form of storage media. A method 
according to the present teachings may also be 
performed by the computational processor that is 
typically part of the VNA. Specifically, it is not 
necessary to know the type of high reflect calibration 
standard 301 used in the calibration process and it is 
not necessary to know the electrical length of the thru 
calibration standard 601. If a phase offset of a 
reflect calibration standard is known, a method 
according to the present teachings is able to measure 
and calculate the electrical length of the thru 
calibration standard 601 rather than relying on 
manufacturer's data. If the electrical length of the 
thru calibration standard is known, a method according 
to the present teachings calculates the phase offset of 
the reflect calibration standard. The calculated 
electrical length or phase offset is more accurate that 
manufacturer supplied data and, therefore, provides for 
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a more accurate calculation of the error artifacts for 
correction of measured data. This is especially 
advantageous when calibrating for on-wafer measurements 
where a repeatable thru calibration standard is 
difficult to obtain and when using a calibration 
standard having characteristics that may have changed 
slightly due to environmental conditions and repeated 
use. In addition, the present teachings recognize that 
calibration using a non-zero through places a 
transmission reference plane at a position different 
from the reflection measurements. A method or 
apparatus according to the present teachings is able to 
calculate the error attributable to use of a non- 
repeatable and non-zero through. The method or 
apparatus according to the present teachings is also 
able to make further corrections to what was considered 
corrected data in the prior art to place the 
transmission measurement reference plane at a same 
point as a position of the reflection measurement 
reference plane. The more accurate through calibration 
standard combined with the additional correction for 
the shifted measurement reference plane provides a more 
accurate measurement of the frequency response 
attributable to only the device under test and not the 
connection environment with which the device under test 
is connected to the measurement equipment. 
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Claims : 

1. A method of measuring a DUT comprising the steps of: 

providing a vector network analyzer having at least 
two measurement ports, 

measuring a reflection characteristic of a high 
reflect calibration standard at each 
measurement port, 

measuring forward and reverse reflection and 
transmission characteristics of a line 
calibration standard, 

measuring forward and reverse reflection and 
transmission characteristics of a source 
terminated through calibration standard, 

measuring forward and reverse reflection and 

transmission characteristics of a locally 
terminated through calibration standard, 

calculating error coefficients for said at least two 
measurement ports based upon results in said 
steps of measuring, 

calculating a shifted electrical length attributable 
to said calibration standards based upon 
results in said steps of measuring, 

measuring a DUT, 

correcting for systematic errors in said step of 
measuring said DUT based upon said error 
coefficients to yield a corrected S-parameter 
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matrix, and 

shifting a reference plane for each element of said 
:rrre"ed S-parameter matrix to coincide with a 
— -asureir.er: t plane. 



, . _ iBrh - d -f measuring as recited in claim - 
" said vectcr network analyzer comprises more 



measurement ports and wherein said steps of 



.„« are repeated for all direct parrs 
measuring are rep«=a ^ 

relP .ent ports and further comprising 



of saia 
the step of 



^ ,,nated through calibration 

..—-«- 1 rn S Ci — — > w 



indirect pairs of measure..*..-. 



3 a method of measuring as recited in claim 2 wherein 
a shifted electrical length between said indirect 
pairs is calculated using load match and source 
match error coefficient terms. 



A me 



thod of measuring as recited in claim 3 wherein 
r 

A „S>1 portn __ C" 
~ ^"M thru nni Lj \2 ihn* nm 

r - - 

A 1.4 _ portm 

wherein S 21 _ thru _™ is equal to S 12 _ thcu _ ro n and an argument 
of both solutions for S 21 _ thr u_™ is fit to a straight 
line, the solution having a y-intercept closest to zero 
being the correct solution and the resulting argument 
of the correct solution being the electrical delay. 
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5. A method of measuring as recited in claim 2 and 
further comprising the step of measuring forward and 
reverse reflection and transmission characteristics 
of a source terminated through calibration standard 
for indirect pairs of said measurement ports. 

6. A method of measuring as recited in claim 2 wherein 
transmission tracking error coefficients are 
calculated using an averaging process for all 
proximal pairs. 

7. A method of measuring as recited in claim 5 wherein 
said step of calculating further comprises 
calculating a different respective shifted 
electrical length for each said direct and indirect 
pair . 

8. A method of measuring as recited in claim 7 wherein 
said shifted electrical length between proximal 
pairs is determined by averaging a shifted 
electrical length between said direct pair and said 
indirect pair having respective proximal pair 
measurement ports in common. 

9. A method of measuring as recited in claim 5 and 
further comprising the step of measuring forward and 
reverse reflection and transmission characteristics 
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of a locally terminated through calibration standard 
for proximal pairs of said measurement ports. 

10. A method of measuring as recited in claim 9 

wherein transmission tracking error coefficients are 
calculated using said locally terminated through 
calibration standard for said proximal pairs. 



11. A method of measuring as recited in claim 2 
wherein said vector network analyzer comprises a 
multiport test set and switch matrix having said 
more than two measurement ports . 

12. A method of measuring as recited in claim 1 and 
further comprising the step of determining a type of 
high reflect calibration standard. 

13. A method of measuring as recited in claim 12 
wherein said step of determining further comprises 
calculating a characteristic of said high reflect 
calibration standard, fitting arguments of two 
possible solutions for said characteristics to a 
straight line, identifying which solution is closest 
to zero phase at DC. 

14. a method of measuring as recited in claim 1 
wherein said ten error coefficients comprise forward 
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and reverse reflection tracking error coefficients, 
wherein said forward and reverse reflection tracking 
error coefficients are determined using a boundary 
condition wherein an argument of reflection tracking 
is zero at DC . 



15. A method of measuring as recited in claim 7 and 

wherein said step of shifting comprises modifying an 
argument of respective S-parameters of said measured 
DUT according to respective ones of said shifted 
electrical lengths. 



16. A method of measuring as recited in claim 15 

wherein said step of shifting a reference plane 
comprises adjusting each said element of said 
corrected S-parameter matrix according to: 

wherein 69 is calculated from said electrical length 
as a function of frequency. 



17. A method of measuring as recited in claim 1 

wherein said step of calculating a shifted 
electrical length comprises calculating a 
characteristic of said high reflect calibration 
standard, fitting an argument of said characteristic 
to a straight line, and using a slope of said 
straight line to calculate a shifted electrical 
length . 



75 



Agilent PDNO 10021239 
inventor: Vahe Adamian 
; Aug 4, 2003 



18. A method of measuring comprising the steps of: 

providing a vector network analyzer having at least 
two measurement ports, 

measuring a reflection characteristic of a high 
reflect calibration standard for each said 
measurement port, 

measuring forward and reverse reflection and 
transmission characteristics of a line 
calibration standard between said at least two 
measurement ports, 

measuring forward and reverse reflection and 
transmission characteristics of a source 
terminated through between said at least two 
measurement ports , 

measuring forward and reverse reflection and 

transmission characteristics of a locally 
terminated through calibration standard between 
said at least two measurement ports, 

calculating forward and reverse reflection tracking 
error coefficients based upon results from said 
steps of measuring using a boundary condition 
wherein an argument of said reflection tracking 
coefficients are zero at DC, 

calculating all remaining forward and reverse error 
coefficients , 

measuring a DUT, 
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cting for systematic errors in said step of 
measuring said DUT based upon said error 
coefficients. 



A method of measuring as recited in claim 18 
wherein said vector network analyzer comprises more 
than two measurement ports and wherein said steps of 
measuring are repeated for direct pairs of said 
measurement ports and further comprising the step of 
measuring a locally terminated through calibration 
standard between indirect pairs of said measurement 
ports . 



20. A method of measuring as recited in claim 18 

and further comprising the steps of calculating an 
electrical length of a shifted reference plane for 
each said at least two measurement ports based upon 
results in said steps of measuring and correcting 
said shifted reference plane to coincide with a DUT 
measurement plane . 



21. A method of measuring as recited in claim 20 

wherein said step of calculating a shifted 
electrical length comprises calculating a 
characteristic of a high reflect calibration 
standard, fitting an argument of said characteristic 
to a straight line and using a slope of said 
straight line to calculate a shifted electrical 
length. 
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22. A method of measuring as recited in claim 19 

wherein said steps of measuring a DUT and correcting 
for systematic errors yields a corrected S-parameter 
matrix and further comprising the steps of 
calculating a respective electrical length of a 
shifted reference plane for each direct pair of 
measurement ports based upon results in said steps 
of measuring and shifting respective elements of 
said corrected S-parameter matrix by respective 
electrical lengths to coincide with a DUT 
measurement plane . 

23. A method of measuring as recited in claim 20 

wherein said step of shifting said shifted reference 
plane comprises adjusting each said element of said 
corrected S-parameter matrix according to: 

S aul =\p\e-^ n) 

wherein 6(6) is calculated from said electrical length 
as a function of frequency. 



A method of measuring as recited in claim 18 
and further comprising the step of determining 
whether said high reflect calibration standard is 
open circuit or a short circuit. 
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25. A method of measuring as recited in claim 24 
wherein said step of determining further comprises 
fitting arguments of two possible solutions for a 
reflection tracking error coefficient to a straight 
line and identifying which solution is closest to 
zero phase at DC. 

26. A method of measuring as recited in claim 19 
wherein a shifted electrical length between said 
indirect pairs is calculated using load match and 
source match error coefficient terms. 



27. A method of measuring as recited in claim 2 6 

wherein 

r 

_ ~~ *J '21 thru nm^ *I2 _thru nm 

*- lA ptmm 

and S 2 i_thru_nm is equal to Si 2 _tnru_mn and an argument of 
both solutions for S 2 i_thru_nm is fit to a straight line, 
the solution having a y-intercept closest to zero being 
the correct solution and the resulting argument of the 
correct solution being used to determine said shifted 
electrical length . 



28. A method of measuring as recited in claim 19 

and further comprising the step of measuring forward 
and reverse reflection and transmission 
characteristics of a source terminated through 
calibration standard for indirect pairs of said 
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measurement ports . 



29. A method of measuring as recited in claim 28 

and further comprising calculating a shifted 
electrical length for each direct and indirect pair 
of measurement ports. 



A method of measuring as recited in claim 29 
wherein transmission tracking error coefficients a 
calculated using an averaging process for all 
proximal pairs of measurement ports. 



31. A method of measuring as recited in claim 28 

wherein said shifted electrical length between 
proximal pairs is determined by averaging a shifted 
electrical length between said direct pair and said 
indirect pair having respective proximal pair 
measurement ports in common. 



32. A method of measuring as recited in claim 28 
and further comprising the step of measuring forward 
and reverse reflection and transmission 
characteristics of a locally terminated through 
calibration standard for proximal pairs of said 
measurement ports . 

33. A method of measuring as recited in claim 32 
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wherein transmission tracking error coefficients are 
calculated using results of measurement taken of 
said locally terminated through calibration standard 
for said proximal pairs. 



34. An apparatus for measuring a DUT comprising: 

a vector network analyzer comprising a signal 

generator and at least two receivers capable of 
connection to at least two measurement ports, 
means for measuring a reflection characteristic of a 
high reflect calibration standard at each 
measurement port, measuring forward and reverse 
reflection and transmission characteristics of 
a line calibration standard, measuring forward 
and reverse reflection and transmission 
characteristics of a source terminated through 
calibration standard, measuring forward and 
reverse reflection and transmission 
characteristics of a locally terminated through 
calibration standard, 

means for calculating ten error coefficients for 

said at least two measurement ports based upon 
results of said means for measuring, 

means for measuring the DUT, 

a processor for calculating a shifted electrical 
length attributable to said calibration 
standards, for correcting systematic errors in 
results from said means for measuring said DUT 
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coefficients to yield a 
matrix, and for shifting 
each element of said 
matrix to coincide with a 



35. An apparatus as recited in claim 34 wherein 

said vector network analyzer comprises more than two 
measurement ports and wherein said means for 
measuring further comprises means for measuring a 
locally terminated through calibration standard for 
indirect pairs of measurement ports. 



36. An apparatus as recited in claim 35 wherein 
said processor calculates a respective shifted 
electrical length between said indirect pairs using 
load match and source match error coefficient terms. 

37. An apparatus as recited in claim 36 wherein 

r 

SA^portn o 

— — 2 1 thru nm ^ 1 2 _ thru nm 

lA_portm 

and S 2 i_thru_nm is equal to S i2 _thru_mn and an argument of 
both solutions for S 2 i_thru_nm is fit to a straight line, 
the solution having a y-intercept closest to zero being 
a correct solution and a resulting argument of said 
correct solution being used to determine said shifted 
electrical delay. 
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38. An apparatus as recited in claim 35 wherein 

said processor calculates transmission tracking 
error coefficients for said proximal pairs using an 
averaging process. 



An apparatus as recited in claim 35 and further 
comprising means for measuring forward and reverse 
reflection and transmission characteristics of a 
source terminated through calibration standard 
between indirect pairs of measurement ports. 



40. An apparatus as recited in claim 39 wherein 
said processor calculates different respective 
shifted electrical lengths for said direct and 
indirect pairs of measurement ports. 

41. An apparatus as recited in claim 40 wherein 
said processor calculates said shifted electrical 
length between proximal pairs by averaging a shifted 
electrical length between said direct pair and said 
indirect pair having respective proximal pair 
measurement ports in common. 



42. An apparatus as recited in claim 39 wherein 

said means for measuring further comprising means 
for measuring forward and reverse reflection and 
transmission characteristics of a locally terminated 
through calibration standard for proximal pairs of 
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said measurement ports. 

An apparatus as recited in claim 42 wherein 
said processor calculates transmission tracking 
error coefficients using results from said means for 
measuring said locally terminated through 
calibration standard for said proximal pairs. 

An apparatus as recited in claim 35 and further 
comprising a multiport test set and switch matrix 
connected to said vector network analyzer for 
providing more than two of said measurement ports. 

An apparatus as recited in claim 34 and further 
comprising means for determining a type of said high 
reflect calibration standard. 

;. An apparatus as recited in claim 45 wherein 

said means for determining further comprises means 
for fitting arguments of two possible solutions for 
a reflection tracking error coefficient to a 
straight line and means for identifying which 
solution is closest to zero phase at DC. 

7. An apparatus as recited in claim 34 wherein 

said ten error coefficients comprise forward and 
reverse reflection tracking error coefficients and 
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said means for calculating determines said forward 
and reverse reflection tracking error coefficients 
using a boundary condition wherein an argument of 
reflection tracking is zero at DC. 



An apparatus as recited in claim 34 wherein 
said measurement ports comprise direct pairs of 
measurement ports and indirect pairs of measurement 
ports and wherein said means for calculating said 
shifted electrical length comprises means for 
calculating respective ones of said shifted 
electrical lengths for each direct and indirect pair 
of measurement ports. 



49. An apparatus as recited in claim 48 and wherein 

said processor for shifting modifies an argument of 
respective S-parameters of said measured DUT 
according to respective ones of said shifted 
electrical lengths. 



50. An apparatus as recited in claim 49 wherein 

said processor for shifting modifies respective 
elements of said corrected S-parameter matrix 
according to: 

S dul =\p\e-'^' n 

wherein 59 is calculated from said electrical length as 
a function of frequency. 
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An apparatus as recited in claim 48 wherein 
said processor calculates said shifted electrical 
length by fitting an argument of a characteristic of 
a high reflect calibration standard to a straight 
line and using a slope of said straight line to 
calculate a shifted electrical length. 



2. An apparatus for measuring comprises: 

a vector network analyzer having at least two 
measurement ports , 

means for measuring a reflection characteristic for 
each said measurement port, measuring forward 
and reverse reflection and transmission 
characteristics for a line calibration standard 
between said at least two measurement ports 
and, measuring forward and reverse reflection 
and transmission characteristics of a source 
terminated through calibration standard between 
said at least two measurement ports, and 
measuring forward and reverse reflection and 
transmission characteristics of a locally 
terminated through calibration standard between 
said at least two measurement ports, 

a processor for calculating forward and reverse 

reflection tracking error coefficients based 
upon a boundary condition wherein an argument 
of said reflection tracking error coefficients 
is zero at DC, and for calculating all 
remaining forward and reverse error 
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coefficients, 
means for measuring a DUT, 

said processor correcting for systematic errors in 
said means for measuring said DUT based upon 
said error coefficients. 



53. An apparatus as recited in claim 52 wherein 
said vector network analyzer comprises more than two 
measurement ports and wherein said means for 
measuring measures all direct pairs of said 
measurement ports and further comprising means for 
measuring forward and reverse reflection and 
transmission characteristics of a locally terminated 
through calibration standard between indirect pairs 
of measurement ports. 

54. An apparatus as recited in claim 52 wherein 
said processor calculates an electrical length of a 
shifted reference plane for each said at least two 
measurement ports and corrects said shifted 
reference plane to coincide with a DUT measurement 
plane . 



55. An apparatus as recited in claim 53 wherein 

said processor calculates respective electrical 
lengths of a shifted reference plane for each said 
direct pair of measurement ports and corrects said 
shifted reference plane for each said direct pair to 
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coincide with a DUT measurement plane. 



An apparatus as recited in claim 55 wherein 
said processor calculates said shifted electrical 
lengths by fitting an argument of a characteristic 
of a high reflect calibration standard to a straight 
line and using a slope of said straight line to 
calculate said shifted electrical length. 



57. An apparatus as recited in claim 55 wherein 

said processor shifts said shifted reference plane 
by adjusting respective elements of a corrected S- 
parameter matrix according to: 

S Jul =\p\e- H ^' n 

wherein 59 is calculated from said respective 
electrical length as a function of frequency. 



58. An apparatus as recited in claim 52, said 

processor also determining a type of said high 
reflect calibration standard. 



59. An apparatus as recited in claim 58 wherein 

said processor determines said type of high reflect 
standard by fitting arguments of two possible 
solutions for a reflection tracking error 
coefficient to a straight line and identifying which 
solution is closest to zero phase at DC. 
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60. An apparatus as recited in claim 53 wherein 
said processor calculates a shifted electrical 
length between said indirect pairs using load match 
and source match error coefficient terms. 

61. An apparatus as recited in claim 60 wherein 

partn ^ r* 

^ 2 1 _ thru nm^\2_ thru _ nm 

l~A _ pnrtm 

and S 21 _thru_™ is equal to S 12 _thru_mn and an argument of 
both solutions for S 2 i_ t hru_nm is fit to a straight line, 
the solution having a y-intercept closest to zero being 
the correct solution and the resulting argument of the 
correct solution being the electrical delay. 

62 . An apparatus as recited in claim 54 and further 

comprising means for measuring forward and reverse 
reflection and transmission characteristics of a 
source terminated through calibration standard 
between indirect pairs of measurement ports. 

63. An apparatus as recited in claim 62 wherein 
said processor calculates a different respective 
shifted electrical length between each direct and 
indirect pair of measurement ports. 

64. An apparatus as recited in claim 63 and wherein 
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said processor for shifting modifies an argument of 
respective S-parameters of said measured DUT 
according to respective ones of said shifted 
electrical lengths. 



An apparatus as recited in claim 64 wherein 
said processor shifts said shifted reference plane 
by adjusting respective elements of a corrected S- 
parameter matrix according to: 

wherein 59 is calculated from said respective 
electrical length as a function of frequency. 



66. An apparatus as recited in claim 63 wherein 
said processor calculates a shifted electrical 
length between proximal pairs by averaging a shifted 
electrical length between said direct pair and said 
indirect pair having respective proximal pair 
measurement ports in common. 

67. An apparatus as recited in claim 63 wherein 
said means for measuring further comprises means for 
measuring forward and reverse reflection and 
transmission characteristics of a locally terminated 
through calibration standard for proximal pairs of 
said measurement ports. 
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An apparatus as recited in claim 67 wherein 
said processor calculates transmission tracking 
error coefficients using results of a measurement of 
said locally terminated through calibration standard 
for said proximal pairs . 



69. An apparatus as recited in claim 63 wherein 
said processor calculates a shifted electrical 
length for proximal pairs of measurement ports from 
an average of shifted electrical lengths for direct 
and indirect pairs having proximal pair measurement 
ports in common. 

70. An apparatus as recited in claim 69 and wherein 
said processor for shifting modifies an argument of 
respective S-parameters of said measured DUT 
according to respective ones of said shifted 
electrical lengths . 



71. A method of measuring a device under test 

comprising the steps of: 

obtaining measured S-parameter characteristics of a 
reflect, a line, and a through calibration 
standard over a range of stimulus frequencies 
for each port of the device under test, 

calculating forward' and reverse reflection tracking 
error coefficients based upon results from said 
measured characteristics using a boundary 
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condition wherein an argument of said 
reflection tracking coefficients are zero at 
DC, 

calculating all remaining forward and reverse error 
coefficients, 

obtaining measured S-parameters characteristics of a 
DUT over said range of frequencies, 

correcting for systematic errors in said measured S- 
parameter characteristics of said DUT based 
upon said error coefficients. 

72. Computer readable media tangibly embodying a 

program of instructions executable by a computer to 
perform a method of determining actual S-parameter 
characteristics of a device under test based upon 
measured S-parameter characteristics of said device 
under test, the method comprising: 

obtaining measured S-parameter characteristics of a 
reflect, line, and through calibration 
standards over a range of stimulus frequencies 
for each port of the device under test, 

calculating forward and reverse reflection tracking 
error coefficients based upon results from said 
measured characteristics using a boundary 
condition wherein an argument of said 
reflection tracking coefficients are zero at 

DC, 

calculating all remaining forward and reverse error 
coefficients, 
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obtaining measured S-parameters characteristics of 
said device under test over said range of 
frequencies , 

correcting for systematic errors in said measured S- 
parameter characteristics of said device under 
test based upon said error coefficients. 
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Abstract 



A method and apparatus comprises providing a vector 
network analyzer having at least two measurement ports 
and presenting a high reflect calibration standard, a 
line calibration standard, a source terminated thru 
calibration standard, and a locally terminated through 
calibration standard at each measurement port or pair of 
measurement ports. Forward and reverse reflection 
tracking is calculated based upon a boundary condition 
wherein an argument of the reflection tracking 
coefficients is zero at DC. The DUT is measured and the 
measurement is corrected for systematic errors based upon 
the error coefficients. 
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